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ION BEAM DIVERGENCE CHARACTERISTICS 


OF TWO-GRID ACCELERATOR-SY-STEMS__ 

Graeme Aston 

* 

Knowledge of Ion beam divergence characteristics for multi -aperture 

accelerator systems is important for both electric space propulsion and 

ground applications of electron-bombardment ion sources. A highly 

collimated beam (minimum off-axis ion velocity components) is generally 

preferred for both applications. The amount of ion beam divergence depends on the 

geometrical grid parameters and operating conditions of the accelerator 

system. Theoretical approaches used to date for predicting the performance 

1-4 

of multiaperture accelerator systems have shown general trends only. 

These theoretical approaches have been limited to the investigation of 
single aperture accelerator systems only. Previous experimental resul-ts 
have shown differences in the ion beam divergence characteristics of single 

5 

and multiaperture accelerator systems. A detailed experimental investi- 
gation was therefore undertaken to determine ion beam divergence over a 

wide range of geometrical and operating parameters of multi aperture accel- 

of 

erat«8» systems. 


Since this section deals only with singly charged positive ions, the 
term "ion” is used with the exclusive meaning "singly charged positive 


ion." 
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Nomenclature 

Jig = Grid separation distance, 
d^ = Screen hole diameter, 
d, = Accelerator hole diameter. 

a 

t^ = Screen grid thickness, 
t, = Accelerator grid thickness. 

a 

A. = Effective acceleration length. 

V-j. = Total accelerating voltage. 

= Net accelerating voltage. 

Vp = Discharge voltage. 

Vg = Sheath potential. 

Vp = Plasma potential. 

R = Net- to- total accelerating voltage ratio, 
a = Ion beam divergence angle, 
fp = Ion beam divergence factor. 

Jg = Beam current per hole, 
q = Ion charge. 

Eq = Free space permittivity. 

= Ion mass. 

= Propellant first ionization potential. 

Ion Beam Formation 

The ion extraction system chosen for this study was the two-grid 
accelerator geometry common to most electron-bombardment ion sources. 
Figure 1(a) depicts a portion of this multiaperture system showing the 
co-axial hole geometry and ion beamlet formation. Figure 2 illustrates 
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D I STANCE 


ACCELERATOR SYSTEM POTENTIAL VARIATION 


FIGURE 2 
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the variation in electrical potential associated with the grid geometry 
of Fig. 1 (a). 

A low density singly ionized plasma is generated within a discharge 
chamber by electron-bombardment of neutral propellant atoms. The dis- 
charge chamber is at a potential a few tens of volts positive of the screen 
grid, which is at a high positive potential (Fig. 2). Around each screen 
hole a stable plasma sheath is formed. This sheath is a space-charge 
limited, emissive surface of positive ions from the discharge plasma. The 

C 

ions reach the sheath with the Bohm velocity” and are accelerated by the 
strong electric field between the positive sheath and negative accelerator 
grid hole (Fig. 2). Depending upon the amount of focusing that occurs, 
most of the ions will pass through the accelerator grid hole and be ejected 
at high velocity into the region beyond (downstream side of the grids). 

The ion, or beam current leaving the accelerator system is space-charge 
limited, so that the accelerator system resembles a vacuum diode. The 

. 7 

perveance for a parallel plate diode can be obtained from Child's law 



The assumption has been made that the area of ion emission is ird|/4. This 
planar expression has been found to agree well with the maximum beam cur- 
rent capacity per hole for a two-grid accelerator system if the effective 
acceleration length is used in place of where 



With this substitution, beam current variations may be expressed in terms 
of normalized perveance per hole 
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NP/H = 




ire 
0 

9 


I 



(3) 


Equation 3 indicates there is an upper limit to the normalized perveance 

per hole for a particular propellant. For mercury, the propellant of 

interest for electric space propulsion, this maximum value is 3.03 x 10"^ 

amp/volt^/^. For argon, the test gas used in this study, this maximum 

value is 6.79 x 10“^ amp/volt^/^. Normalized perveance per hole 
/ 3 / 2 \ 2 

rB/Vf extraction performance 

of each grid set investigated. 

Ion focusing effects occur in two regions of the accelerator system; 
between the grids and just downstream of the grids. Ions between the grids 
tend to follow the electric field lines passing from the screen hole plasma 
sheath through the accelerator hole. This tendency to follow electric 
field lines is most evident near the plasma sheath, where the ion velocity 
is small. The position and shape of the screen hole plasma sheath is crucial 
to this phase of the ion focusing process because it determines the local 
accelerating electric field distribution. Both screen hole plasma sheath 
position and shape depend on the ion current being extracted and the accel- 
erator system geometry. In normal operation the shape of the plasma sheath 
is similar to that shown in Fig. 1(a). As plasma density increases the 
sheath moves further into the screen hole to satisfy the requirement of 
space-charge limited ion current emission from the discharge chamber plasma. 
The sheath shape becomes flatter and direct ion impingement upon the accel- 
erator grid begins to occur (Fig. 1(b)). Further plasma density increases 
result in further ion focusing degradation. In practice this high impinge- 
ment limits the ion current which can be extracted from an accelerator 


system. 
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As the focused ions leave the accelerator system, electrons are in- 
jected into the beam to maintain a neutral charge efflux. These neutral- 
ization electrons spread rapidly throughout the beam producing a neutral 
plasma. Because of the negative accelerator grid potential, neutralization 
electrons which approach the accelerator holes are repulsed by electro- 
static forces. The result is a neutralization, or equipotential , surface 
created slightly downstream from the accelerator grid which is at approxi- 
mately facility ground potential (Fig. 1(a)). Between the accelerator 
holes and the neutralization surface the ions are decelerating through a 
region of retarding electric field (Fig. 2). Although the deceleration 
incurred is less than the acceleration imparted by the grid set, it is still 
sufficient to defocus the ion trajectories significantly. However, once 
these ions reach the neutralization surface, electrostatic shielding by the 
beam plasma nullifies the accelerator grid effect and little additional 
change in ion trajectories occurs. Ion defocusing decreases as the net-to- 
total accelerating voltage ratio R, (Fig. 2) increases so it would appear 
desirable to have R as near unity as possible. In practice however, a 
value of R = 1.0 cannot be reached because neutralization electrons will 
backstream through the accelerator holes when the negative accelerator 
grid potential is too small. 

Both regions of ion focusing mentioned above were examined for their 
contribution to the overall ion beam divergence. Ion focusing between the 
grids was varied by changing the ion beam current, the accelerator system 
geometry, and the total accelerating and discharge voltages. Ion defocusing 
was varied by changing the net-to-total accelerating voltage ratio. 
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Apparatus and Procedure 

For this study, a mildly divergent magnetic field 8-cm diameter 
electron-bombardment ion source was constructed. Tungsten wire filaments 
were used as both the main and neutralizer cathode emitters. The magnetic 
field was derived from a long solenoidal winding extending the length of 
the discharge chamber, with additional coil windings positioned at the 
chamber's upstream end. The coil geometry was such that the field at the 
downstream end was sixty percent that of the upstream end. A cylindrical 
anode was employed and non-magnetic stainless steel construction used 
throughout the source. Further design details and operating characteris- 
tics for the source are found in Reference [8]. Mercury is the propellant 
of principle interest for electric space propulsion, but because of the in- 
ability to change ion source operating conditions rapidly with mercury it 
was not used for this study. Argon was chosen as propellant because of 
its availability and the ease and speed with which flow rates could be set 
and altered. 

The screen and accelerator grids were made from thin flat sheet 
graphite. The grid aperture pattern used comprised a nineteen hole 
hexagonal array with a center- to- center hole spacing of 2.54 mm. It was 
felt this number of holes adequately modelled the adjacent screen hole 
interactions found in full-size grid systems. Variable grid separation 
was accomplished by using various numbers of thin mica washers (0.254 mm 
thick). The grid sets were fastened together with four insulated stainless 
steel, bolts. Alignment of the screen and accelerator grids was accomplished 
by hand beneath a large illuminated magnifying glass; this straightfor- 
ward technique was found to produce good data reproducibility. The assem- 
bled accelerator system was placed on a masked down discharge chamber. 
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which insured ion extraction from a near uniform plasma. Calculations 
predicted the greatest spatial variation in discharge chamber plasma 
density across the 1 cm^ hexagonal hole array would be less than 5%. 

Grid separation could be checked with feeler gauges when cold, but while 
the source was operating and the grids were hot they could only be ex- 
amined visually. While some grid warpage could occur during operation, 
the magnitude of this warpage should be small because the grids were made 
of carbon and the greatest distance between grid supports was less than 
2-cm. All source operation was conducted in a 30-cm pyrex bell jar. 

Average bell jar pressure was 1 x 10“^^ torr. 

Beam divergence measurements were accomplished using a movable probe 
rake containing twenty individual Faraday ion current sensors described 
in Reference [8]. This probe rake was positioned normal to the beam axis 
and 16.7 cm downstream of the accelerator system. Figure 3 shows the ion 
source and probe rake in position. The probes were screened and biased 
to reflect neutralization electrons and"low energy charge exchange ions 
present in the beam plasma. After stable source operation was achieved, 
the probe rake was centered along the beam axis. Beam current was varied 
by adjusting the refractory cathode current and hence electron emission 
and discharge chamber plasma density. Each grid set geometry investigated 
was operated over a range of increasing emission levels to the maximum 
beam current obtainable from that grid set. The approach to this maximum 
was usually characterized by a rapid increase in accelerator grid impinge- 
ment current and a negligible beam current increase. Ion beam profiles 
were obtained at selected beam currents leading up to this maximum. Each 
profile consisted of twenty ion current density measurements, corresponding 
to the ion current received by each Faraday probe as a function of its 
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METHOD OF OBTAINING ION BEAM PROFILES 


FIGURE 3 
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position normal to the beam axis. Figure 3 shows a typical ion beam 
profile. From these profiles ion beam divergence could be quantified 
by calculating the ion beam divergence angle a, (Fig. 3). The ion beam 
divergence angle was defined as the truncated cone angle enclosing 95% 
of the total integrated beam current. Another less obvious parameter 
for determining ion beam divergence is the ion beam divergence factor f^, 
which is of considerable interest for ion thruster applications of 
electron-bombardment ion sources. The divergence factor is the ratio of 
the net axial thrust produced by the divergent ion beam to -the thrust that 
would be produced if the ion beam were perfectly collimated. Both ion 
beam divergence angle and ion beam divergence factor were the calculated’ 
parameters (based on experimental measurements) of primary interest to 
this study. The following equations were used in their determination 


R - R 

a = Arctan ( j) 


f 


D 


77 

S J. R. Cos3. 
i=l 1 T 1 


E J. R. 


(4) 


(5) 


where 

R^ = radius normal to the beam axis defining a cone enclosing 
95% of the total integrated beam current, 

Rg = ion beam radius at accelerator grid , 

= probe current for the i^^ probe (an interpolation routine 
was used to increase the number of probe data points from 
20 to 77) ^ 

*fch 

R^. = distance from ion beam axis to i^" probe , 
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3^. f angle formed by intersection of a line extending from the 
center of the, ion source grid set to the i^*^ probe and the 
ion beam axis. 

Results and Discussion 

The effect of grid geometry variations on ion beam divergence was 
determined by independently varying the grid separation distance 
accelerator hole diameter d , screen grid thickness t. and accelerator 

a 5 

grid thickness t,. To non-dimensional ize the geometrical grid parameters 

a 

they were divided by the screen hole diameter d^, which remained at 
2.06 mm throughout this study. The precision of the data obtained was 
verified by repeatedly testing various grid set geometries at different 
times with the grids reassembled and realigned for each test. The results 
of these tests showed a maximum variation in a and f^ of t 0.5® and 
t 0.001 respectively. The absolute accuracy of the data is thought to be 
better than t 1.0® for a, and better than t-0.002 for f^. A tabular listing of 
the experimental results presented graphically in this section can be found 
in Appendix I. 

Figure 4 shows the effect of propellant flow rate variations on ion 
beam divergence for a grid geometry, the total accelerating voltage, dis- 
charge voltage and net-to-total accelerating voltage ratio typical of 

910 " 1112 

present mercury ion thrusters ’ and ground-based argon ion sources. ’ 

An almost five-fold variation in flow rate is seen to alter the grid 

systems divergence characteristics to a very minor extent. The small, 

consistant beam divergence angle decrease with increasing propellant flow 

1 3 

rate shown in Figure 4 has been reported previously by Altenburg, et.al. 

The present data verifies those observations. For very small flow rates 



- 13 - 



EFFECT OF ARGON FLOW RATE ON ION BEAM DIVERGENCE 


FIGURE 4 
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the discharge became unstable and extinguished before any accelerator 
system limitation was reached (circles in Fig. 4). For the remainder of 
the tests reported herein the ion source was operated at the minimum flow 
rate which would give stable source operation up to the maximum beam cur- 
rent obtainable from the grid set being tested. This resulted in a pro- 
pellant flow rate which varied from 0.075 to 0.125 Standard Cubic Centi- 
meters per Minute (S.C.C.M.) over the range of grid set geometries and 
operating conditions investigated. Figure 4 shows these flow rate varia- 
tions should produce a negligible effect on the ion beam divergence 
characteristics of a particular grid set geometry. 

Figure 4 shows also that there is an optimum perveance, corresponding 
to an optimum screen hole sheath position and shape, which gives the lowest 
ion beam divergence angle. This minimum divergence angle defines the accel- 
erator system operating condition most desirable from an ion focusing view- 
point. It is interesting to note that grid set operation at the minimum 
beam divergence angle is also close to the maximum beam current at which 
one would prefer to operate. Increasing the beam current beyond this point, 
thereby allowing the screen hole sheath to become flatter in shape (Fig, 1(b))., 
results in an increasing fraction of accelerated ions impinging directly 
upon the accelerator grid'. The assocfated increase in ion sputtering 
damage seriously degrades accelerator grid lifetime. 


Figures 5 through 8 show the beam divergence variations for changes to 

the grid separation ratio , and net-to-total accelerating voltage ratio R, 

°s d 

For these data the accelerator hole diameter ratio -A , screen grid thickness 

‘s = 

ratio ^ , and accelerator grid thickness ratio , were held constant at 

°s °s 

values of 0.642, 0.185 and 0.370 respectively. Qualitatively all curve 


shapes are similar. This similarity in curve shape and the fact that 
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FIGURE 5 
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★ 

variations in the maximum normalized perveance per hole are small indicates 
the usefulness of normalized perveance per hole as a correlating parameter. 
Increasing the grid separation ratio for a constant net-to-total accelerating 
voltage ratio decreases ion beam divergence angles, with corresponding ion 
beam divergence factor increases. This is presumably a result of the more 
planar equipotential surfaces associated with the larger grid separation 
distance. Similarly, increasing the net-to-total accelerating voltage ratio 
for a constant grid separation ratio decreases overall ion beam divergence. 
This effect is greatest for the smallest grid separation ratio tested and 
•decreases with. grid' separation ratio increases. The uniform manner in which 
net-to-total accelerating voltage ratio variations raise or lower the ion 
beam divergence angle curves of Figs. 5 through 8, suggests these ion de- 
focus ing variations do not affect the ion focusing processes characteristic 
of each grid separation ratio. Figures 5 through 8 show that large ion beam 
divergence variations are possible with two-grid accelerator systems. A 

change in the minimum ion beam divergence angle from 25.1° to 4.9° is 

£ ■ £ 

observed' when going from values of = 0.247 and R = 0.5 to ^ = 1.000 
and R = 0.9. 

Figures 9 through 11 show the effect of accelerator hole diameter 
ratio, screen grid thickness ratio and accelerator grid thickness ratio 
variations on ion beam divergence. Beam divergence angle changes are 
observed to be very slight for large variations in these parameters. 

However, decreasing the accelerator hole diameter ratio substantially de- 
creases the maximum normalized perveance per hole, or beam current. 

It is important to realize that constant normalized perveance per hole 
does not mean constant beam current. If the total accelerating voltage 
Vj, is held constant, the beam current will increase as the grid 
separation ratio is decreased (Equation 3). 







EFFECT OF SCREEN GRID THICKNESS RATIO 
ON ARGON ION BEAM DIVERGENCE 

FIGURE TO 
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EFFECT OF ACCELERATOR GRID THICKNESS RATIO 
ON ARGON ION BEAM DIVERGENCE 


FIGURE n 
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obtainable from the grid set (Fig. 9). The larger accelerator hole 
apparently allows the screen hole plasma sheath to penetrate further into 
the screen hole, with reduced ion focusing (Fig, 1(b)), before an ex- 
cessive accelerator grid ion impingement current is reached. Larger 
accelerator holes result therefore In increased beam currents. Of particular 
interest in Fig. 9 are the results obtained with the smallest accelerator 
hole diameter ratio tested. Here the beam divergence angle goes through a 
maximum with increasing beam current. The accelerator grid impingement cur- 
rent observed with this grid set was large at the smallest beam current, 
went through a minimum as beam current was increased and then rose sharply 

as the maximum beam current was approached. These characteristics suggest 

14 

ion trajectory cross-over may be occurring at low beam current levels. 

If this is occurring the apparent decrease in divergence angle at very low 
normalized peryeances could be caused by an increase in the number of ions 
leaving the screen hole sheath edge that are impinging directly upon the 
accelerator grid. In this case the accelerator grid would be acting as a 
beamlet mask for severe off-axis velocity ions at very small accelerator 
hole diameter ratios. 

The different grid set geometries presented in Figs. 9 through 11 were 
investigated for the effects of net-to-total accelerating voltage ratio 
variations also. No significant changes in the shape of the ion beam di- 
vergence angle curves, found characteristic to each geometrical grid parameter, 
were apparent. Rather these curves were simply raised or lowered by a 
uniform amount, predictable from Figs. 5 through 8, as the net-to-total 
accelerating voltage ratio was varied. 

In summary. Figs. 5 through 11 show that the geometrical grid parameters 
and net-to-total accelerating voltage ratio act fairly independently of one 
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another in the effect each has on the ion beam divergence 'characteristics 
of two-grid accelerator systems. These effects can be summarized as 
fol 1 ows ; 

(i) Ion beam divergence depends most strongly on grid separation 

ratio ) and net-to-total accelerating voltage ratio (R) 

9 5 

changes. 

(ii) The maximum normalized perveance per hole (or beam current) 

obtainable from a particular grid set depends most strongly on 
the accelerator hole diameter ratio (d /d ). 

a S 

(iii) Screen grid thickness ratio (t^/d^) and accelerator grid thickness 

ratio (t /d ) variations result in consistent, but negligible, 
a s 

changes in ion beam divergence and the maximum normalized 
perveance per hole. 

While these trends were obtained for a total accelerating voltage and dis- 
charge voltage of 1100 and 40 volts, respectively, the same trends have 
been reported previously for a total accele“fating voltage and discharge 
vol tages of 600 and 40 volts respectively. This indicates their general nature. 

Figures 12 through 14 show the effect of total accelerating voltage 
and discharge voltage variations on ion beam divergence. For these data 

the net-to-total accelerating voltage ratio was kept at a constant value 

n 

of 0.7 while the grid geometry parameters remained at the values = 0.741, 
d t t °s 

= 0.642, = 0.185 and = 0.370. Increasing the total accelerating 

s . s s 

voltage (for a constant discharge voltage) or decreasing discharge voltage 
(for a constant total accelerating voltage) shifts each beam divergence 
curve to higher normalized perveance per hole values. For example, the 
normalized perveance per hole at which a beam divergence angle of 15° occurs 
is larger at = 1500 'volts than Vj = 600 volts (for constant discharge 
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EFFECT OF DISCHARGE VOLTAGE ON ARGON ION BEAM DIVERGENCE 
(Total accelerating voltage - 600 volts) 


FIGURE 12 
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voltage) and is larger also at = 32 volts than Vp = 50 volts (for 
constant total accelerating voltage). 

The changes in the minimum ion beam divergence angle are small, but 
there is a consistant trend towards increased beam divergence with in- 
creasing total accelerating voltage and decreasing discharge voltage. 

This effect is most pronounced for V.j. = 600 volts (Fig. 12) where a 3° 
average beam divergence angle increase is observed for a discharge voltage 
decrease from 50 to 32 volts. For higher total accelerating voltages of 
1100 and 1500 volts (Figs. 13 and 14) the angular spread caused by dis- 
charge voltage variations is less. The si ight reduction in maximum normalized 
perveance per hole as Vj^ is decreased from 40 to 32 volts for total 
accelerating voltages of 1100 and 1500 volts (Figs. 13 and 14), is thought 
to be caused by a discharge power supply limitation when operating at 
these extreme conditions. 

Variations in beam divergence with discharge voltage are not con- 
sidered to be due to changes in doubly-charged ion densities because these 
densities should be negligible in this 8 cm argon thruster at discharge 
voltages below 50 v.^^ Rather discharge and total accelerating voltage 
variations are believed to alter the sheath shape and position. The manner 
in which these changes occur to produce the results shown in Figs. 12 
through 14 however is still open to conjecture. It can be said though that 
these changes have less effect on the beam divergence characteristics of a 
particular grid geometry as the discharge- to- total accelerating voltage 
ratio is decreased. Figure 15 combines beam divergence data from four 
different total accelerating voltages using a constant 40 volt discharge. 

The geometrical grid parameters and net-to-total accelerating voltage ratio 
remained at the values indicated. Figure 15 shows the beam generally be- 
comes less divergent as the discharge- to- total accelerating voltage ratio 
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EFFECT OF TOTAL ACCELERATING VOLTAGE ON ARGON ION BEAM DIVERGENCE 
(Discharge voltage - 40 volts) 


FIGURE 15 
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is decreased. Increasing total accelerating voltage above about 1500 
volts alters ion beam divergence by a negligible amount only. Theoretical 

4 

results obtained by Kaufman are shown also. These theoretical results 
were based on a zero discharge voltage and this is believed a major reason 
for their inability to model the experimental results adequately. Figure 15 
indicates closer agreement between experiment and theory as the experi- 
mental discharge-to-total accelerating voltage ratio approaches zero. 

The data contained in Figs. 12 through 15 show that neither total 
accelerating voltage nor discharge voltage alone adequately describe the 
beam divergence variations. From the observed opposite trends that in- 
creases in Vt and produce, together with dimensional considerations, 

the ratio of these two parameters ^ , is believed to be the most important 

T 

quantity to describe the overall effect. Other grid set geometries 

z 

(different etc.) were also operated to determine the effect of discharge 
s 

and total accelerating voltage variations on divergence characteristics. 

No differences in the voltage induced trends depicted in Figs. 12 through 
15 were observed. The only differences apparent were related directly to 
the grid geometry changes. For instance, a uniform increase in ion beam 
divergence occurred after decreasing the grid separation ratio from 0.741 
to 0.494, but the amount was predictable from the results shown in Figs. 6 
and 7. 

Correlated Results 

In the previous section it was shown that ion beam divergence, from 
two-gridj multi -aperture accelerator systems, depends primarily on the grid 
separation ratio , net- to- total accelerating voltage ratio R, and dis- 
charge vol tage-to-total accelerating voltage ratio ♦ These three 
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parameters have been used to correlate the ion beam divergence angle 

variation for all grid set geometries and accelerator system operating 

conditions examined during this study. Figure 16(a) shows the result 

of this correlation. Here the simple Child's Law normalized perveance 

per hole expression has been modified to include the effect of discharge- 

to- total accelerating voltage ratio variations. The ion beam divergence 

angle correlation is seen to lie within a band t H° (actual degrees) of 

the indicated line of best fit. The parameter NP/Hq . occurring in the 

min 

beam divergence angle correlation refers to the normalized perveance per 
hole at the minimum beam divergence angle. It is plotted as a function 
of the discharge vol tage-to-total accelerating voltage ratio ^ , and 
shown in Fig. 16(b). It was necessary to introduce normalized perveance 
per hole as a parameter in -the ion beam divergence angle correlation be- 
cause of the decreased .effect net-to-total accelerating voltage ratio 
variations have on ion beam divergence angle at low normalized perveance 
per hole values (Figs. 5 through 8.). Figures 16(a) and 16(b) allow ion 
beam divergence angles to be determined for a wide variety of grid geometry 
and accelerator system operating conditions. However, as mentioned prev- 
iously, there is an absolute beam current limit, related directly to the 
accelerator hole diameter ratio, for each grid set geometry. This limit 
is shown graphically in Fig. 16(c). Figure 16(c.) may be interpreted as 
describing the beamlet diameter ratio at the accelerator hole entrance as 
a function of the correlated normalized perveance per hole, because the 


beamlet diameter is approximately equal to the accelerator hole diameter 

at the maximum normalized perveance per hole (or beam current) condition. 

Figures 16(a), 16(b), and 16(c) should be used in the following way. 

^e 

Firstjthe values of Jg, , V^, V.^. and R must be known. The value of 
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NP/Hco^^ is then calculated and used to obtain the value of ■from 

Fig, 16(a). The value of NP/H„ is found from Fig. 16(b) and used to 

“min 

calculate the actual beam divergence angle from the previously found 

value of Finally, Fig. 16(c) should be used to determine the 

minimum value of d to attain the desired beam current. To avoid excessive 

a 

accelerator grid ion impingement this minimum d. value should be increased 

a 

by about twenty percent. If accelerator system operation over a range of 
beam currents is proposed. Fig. 16(c) should be used to determine the value 
of d^ (by the method outlined above) necessary to attain the maximum beam 

a 

current desired. The procedure outlined previously will then establish the 
range of divergence angles to be expected for the beam current operating 
range. 

Figure 17 shows the minimum ion beam divergence angle that can be 
obtained as a function of grid separation ratio and net-to-total accel- 
erating voltage ratio. A straight line was chosen as a conservative fit 
to the extrapolated regions which are~considered to be good for a first 
approximation of the minimum beam divergence angle when operating at these 
extreme conditions. The data shown in Fig. 17 is most accurate for the 
discharge voltage and total accelerating voltage indicated. However, prev- 
ious results (Fig. 12 through 15) have shown that varying and either 

side of 40 and 1100 volts respectively, produced a variation of t only 
in the minimum beam divergence angles given in Fig. 17. Hence Fig. 17 may 
be used to obtain fairly representative minimum beam divergence angle 
trends for a wide range of discharge- to- total accelerating voltage ratios. 

Screen Hole Size and Propellant Variations 

Figures 18 and 19 show the results of varying the two parameters held 
constant throughout this study, screen hole size and propellant gas. 
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A constant grid aperture array size, similar non-dimensionalized grid 
geometry parameters and the same argon propellant mass flow rate were 
used for each screen hole size. Similar discharge chamber neutral number 
densities with the same grid system geometry were used for the different 
propellants (an exception was neon which required twice the neutral 
number density of the other gases, a 46 volt discharge and a discharge 
chamber magnetic field reduction before a stable discharge plasma could 
be maintained). Normalized perveance per hole for the different propellants 
shown in Fig. 19 was corrected to correspond to operation with argon. 

This was done by multiplying each calculated normalized perveance per hole 

value (Equation 3} by the square root of the ratio of the propellant atomic 

1 

weight to the atomic weight of argon, (at. wt. prop. /at. wt. argon)^. The 

reduction in maximum normalized perveance per hole with decreasing screen 

hole size has been observed previously with two-grid electron-bombardment 
14 

ion thrusters. The present data verifies those observations but does 
not shed new light onto the cause of'the effect. 

Figure 19 shows the parameter, normalized perveance per hole, corre- 
lates successfully the ion beam divergence variation over a range of pro- 
pellant atomic weight and ionization potential. Figures 18 and 19 show 
that screen hole size and propellant gas changes do not alter substantially 
the beam divergence variation found characteristic of a particular grid set. 
Therefore, the data and correlation parameters presented in this report, 
while being most accurate for argon, are considered to hold generally for 
other propellants and screen hole sizes also. For instance, if knowledge 
of the ion beam divergence characteristics of a particular grid geometry 
operating with mercury rather than argon propellant are desired, then the 
following procedure should be used. From Equation (3) the normalized 
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perveance per hole is calculated for the grid set operating condition 
of interest. This calculated normalized perveance per hole value is 
multiplied by the square root of the ratio of the atomic weight of 

X 

mercury to the atomic weight of argon, (at. wt. mercury/at. wt. argon)^ 

= 2.241. Now, having corrected the original value of normalized perveance 
per hole for mercury operation to a normalized perveance per hole value 
for use with argon, the figures and correlation parameters presented in 
this report may be used. 

Conclusions 

A comprehensive experimental investigation of two-grid accelerator 
systems has been conducted. A wide range of geometrical grid parameters 
and grid set operating conditions were investigated for their effect on 
ion beam divergence. Ion beam divergence was found to depend most strongly 
on grid separation ratio, net-to-total accelerating voltage ratio and dis- 
charge- to-total accelerating voltage ratio variations. The maximum 
normalized perveance per hole (and hence beam current) obtainable from a 
particular grid set, at a constant discharge-to-total accelerating voltage 
ratio, was found to depend most strongly on the accelerator hole diameter 
ratio. Variations to the other geometrical grid parameters have a 
negligible effect on ion beam divergence. Graphical results were obtained 
which permit performance predictions in new accelerator system designs. A 
general divergence angle correlation was developed to permit approximate 
divergence estimates of present designs. All results contained in this 
paper are considered of sufficient reliability to allow interpolation 
between specific experimental results and the use of limited range extra- 
polations. Although argon was the main test gas used in this study, the 
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results are applicable to other propellants when the correction to 
normalized perveance per hole has been made (see pages 39-40). 



ION BEAM DIVERGENCE CHARACTERISTICS OF THREE-GRID 
ACCELERATOR SYSTEMS 

Graeme Aston 

The addition of a third (decelerator) grid to the conventional 
two-grid accelerator system is known to reduce ion beam divergence.^ ^ 

The decelerator electrode is maintained at, or close to, facility ground 
potential and positioned downstream of the accelerator grid. This fixes 
the neutralization surface at the decelerator grid potential, resulting 
in flatter equipotential surfaces between the neutralization plane and 
accelerator hole than with a two-grid accelerator system. Flatter 
equipotential surfaces in the ion deceleration region cause less off-axis 
deflection of the emerging ions. This results in decreased ion beam 
divergence over that obtained by using a two-grid accelerator system. 

The purpose of this study has been to determine those geometrical 
and operating conditions which affect moststrongly the ion beam divergence 
of three-grid accelerator systems. The two-grid accelerator system studies 
in the preceding section have shown that screen and accelerator grid thick- 
ness and hole diameter changes affect ion trajectories in the ion decelera- 
tion region only slightly. Hence, decelerator grid location and geometry 
and accelerator system operating conditions have been the parameters in- 
vestigated most thoroughly in this report. 

Parameter Definitions 

Figure 20 shows a portion of a multi aperture three-grid accelerator 
system indicating the geometrical parameters used. In the data presented 
hereafter these parameters will have the following values unless specified 
otherwise: 
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(i) Screen hole diameter, = 2.06 mm, 

(ii) Screen grid thickness ratio, t^/d^ = 0.185, 

(iii) Accelerator grid thickness ratio, t /d = 0.370, 

0. S 

(iv) Decelerator grid thickness ratio, tyd^ = 0.370, 

(v) Accelerator hole diameter ratio, d /d = 0.642, 

a S 

(vi) Decelerator hole diameter ratio, d^/d^ = 0.827, 

(vii) Total accelerating voltage, = 1100 volts, 

(viii) Discharge voltage, = 40 volts, 

Perveance is based on the argon propellant ion mass because argon was used 
in the tests. Results are considered correct for other propellant atoms 
if perveance values are multiplied by the square root of the ratio of argon 
atomic mass number (40) to that of the propellant to be used. To modify 
the data of this section for use with mercury propellant for example 
perveance values, given on the axis of the graphs would be multiplied by 
/40/200 = 0.447. Perveance and other parameters used in this section are 
calculated as discussed in the preceding section of this report. 

Apparatus and Procedure 

The apparatus and procedure used for this study were identical to 
that which has been discussed in the preceding section of this report 
except for the addition of the decelerator grid. A twelve volt negative 
bias was placed on this grid to prevent neutralization electrons from 
reaching its surface. Both accelerator and decelerator grid impingement 
currents were monitored for each grid set operating condition. The pre- 
cision with which the ion beam divergence data were obtained is t 0.5° 
for the ion beam divergence angle a, and t 0.001 for the ion beam diver- 
gence factor fp. The absolute accuracy of the data is thought to be better 
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than t 1° for a, and better than t 0.002 for fp. 


Results and Discussion 

Figures 21 through 24 show the beam divergence characteristics of a 
high peryeance grid set (close screen-to-accelerator grid separation ratio 

) as the accelerator-to-decelerator grid separation ratio and net-to- 
s • , °s 

total accelerating voltage ratio R are varied over the range of interest. 

Either increasing the net-to-total accelerating voltage ratio or decreasing 

the. accelerator-to-decelerator grid separation ratio is observed to decrease 

ion beam divergence. The ion beam divergence becomes less sensitive to 

these parameters as the net-to-total accelerating voltage ratio approaches 

unity. Figures 21 and 22 show that reducing the accelerator-to-decelerator 

grid separation ratio by a factor of two from 0,247 to 0.123 does little to 

lower the minimum beam divergence angle. This indicates that a limiting 

beam divergence angle minimum is approached which seems unaffected by further 

accelerator-to-decelerator grid separation ratio reductions beyond a 

°s 

value of 0.2 (for this particular grid set). 

The beam divergence angle reductions at low normalized perveance seen 
in Figs. 21 through 24, result from direct ion impingement upon the decel era- 
tor grid — a condition that results in decel grid erosion and probably makes 
operation at low perveance impractical. This reduction in divergence angle 
at low perveance appears to occur because the beamlet ions impinging upon 
the decelerator grid at low beam current levels are those which have the 
greatest off-axis velocity components. By intercepting these ions the decelera- 
tor ’grid acts -as a beamlet mask which reduces the beam divergence angle at 
low beam current -levels.- Evidence of the high decelerator impingement at 
low perveance is shown in Figure 25. Here accelerator and decelerator 
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impingement currents, expressed as a percentage of beam current are 
plotted against normalized perveance, for the extreme cases of accelerator- 
to-decelerator separation presented in Fig. 21. At low beam currents (low 
normalized perveance) the decelerator grid impingement current is high. As 
beam current increases the decelerator grid impingement current drops 
quickly to a minimum value and then- increases as the maximum beam current 
for that particular grid set is approached. By contrast, the accelerator 
grid- impingement current increases slowly with beam current until the maxi- 
mum beam current is approached, where it too increases sharply. Decelerator 
impingement is observed to increase sharply with an increase in accelerator- 
to-decelerator grid separation (^) (Figure 25a). Figure 25b shows the 
decelerator grid impingement current variation as a function of the net-to- 
total accelerating voltage ratio. Variations to the net-to-total accelerating 
voltage ratio, for a given grid set geometry, are seen to affect the decel- 
erator grid impingement current less strongly than accelerator- to-decelerator 
grid separation ratio changes. 

Figure 26 is a plot of minimum beam divergence angle obtained from 
Figures 21 to 24 against net-to-total accelerating voltage ratio, with 
accelerator- to-decelerator grid separation ratio as a parameter. A two-grid 
accelerator system divergence angle curve is also shown for comparison. 
Significant beam divergence angle reductions for three-grid accelerator 
systems over two-grid accelerator systems are apparent, even for a rela- 
tively high net-to-total accelerating voltage ratio of 0.7. Extrapolating 
the two and three-grid minimum divergence angle curves of Fig. 26 shows 

t 

that they'cross at a net-to-total accelerating voltage ratio of about unity. 
Figure 27 compares minimum' beam divergence angle results obtained-experi- 
mentally to .those 'Obtained theoretically by Kaufman. Qualitative agreement 
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between theory and experiment is seen to be good, but the experimental 
results do show a more substantial reduction in divergence angle associ- 
ated with the substitution of three for two-grid accelerator systems than 
the theory predicts. 

Perhaps the most striking difference between the results of Figs. 21 
through 24 and two-grid data contained in the preceding section lies in 
the very sharp approach to the minimum divergence angle for three-grid 
accelerator systems. Typically, the operation of three-grid systems at, 
or near, the minimum beam divergence angle is restricted to a very narrow 
normalized perveance (or beam current) range. Unfortunately, this problem 
is most severe for low net- to- total accelerating voltage ratios, where 
three-grid accelerator systems have their greatest focusing advantage and 
hence greatest potential use. Although, as was previously mentioned, only 
negligible focusing gains are realized by reducing the accelerator-to- 
decelerator grid separation ratio ^^d^ from 0.247 to 0.123, Figs. 21 and 
22 do show a significant increase in the mininfum divergence angle operating 
range is achieved by doing this. Successful three-grid operation would 
require accurate placement of the accelerator system operating condition 
at the minimum beam divergence angle operating point. Failure to do this 
would result in accelerator system operation on the steeply varying portion 
of the divergence angle curve. Here, small beam current fluctuations would 
result in large beam divergence angle variations; a condition which could 
result in low integrated thrust unless the ion source were beam current 
regulated. 

Figure 28 shows the results of preliminary tests designed to determine 
the best decelerator hole size and accelerator-to-decelerator grid separa- 
tion ratio for a given screen and accelerator grid geometry. For, these 
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tests a decelerator hole diameter ratio of d,/d = 0.827, and accelerator- 

d s 

to-decelerator grid separation ratio of = 0.247, gave the lowest ion 

beam divergence. Figure 29 shows that significant decelerator grid impinge- 
ment occurs as the decelerator hole size approaches the accelerator hole 
size. For equal accelerator and decelerator hole sizes the decelerator 
grid impingement current is unacceptably high. These initial results in- 
dicate ’that low beam divergence and low decelerator grid impingement 
currents for three-grid accelerator systems require that the decelerator 
hole diameter be about twenty percent larger than the accelerator hole 
diameter (and the decelerator grid be as close to the accelerator grid as 
possible, as was pointed out earlier). Further work is being done to 
validate these accelerator system design criteria for a larger range of 
grid geometries and operating conditions than those shown in Fig. 28. 

It should be noted that small negative decelerator grid impingement 
currents were observed whenever the beamlet diameter was substantially 
less than the decelerator hole diameter. This occurs for large decelerator 
hole diameter ratios, small accelerator- to-decelerator separation ratios 
and large net-to-total accelerating voltage ratios. These negative decel- 
erator grid impingement currents are thought to be a result of secondary 
electrons produced by direct ion and charge exchange ion impingement upon 
the accelerator grid. The secondary electron currents are always present 
but they are measureable only when the ion impingement upon the decelerator 
grid approaches zero. These negative decelerator grid impingement currents 
were small, being at most 0.2% of the beam current. 

Figure 30 shows the beam divergence characteristics of a low perveance, 

good focusing grid set (large screen- to-accelerator grid separation ratio 

s, / 4 ), as the accelerator- to-decelerator grid separation ratio Aj/d , is 
g s US 
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Vj= 1 100 volts 
Vn= 40 volts 
R = 0.3 





EFFECT OF DECEL GRID HOLE DIAMETER ON 
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FIGURE 29 
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varied for a constant net- to- total accelerating voltage ratio of 0.3. 
Comparing Figs. 21 and 30 shows some striking differences. The most 
notable being that the minimum beam divergence angle in Fig. 30 occurs 
not for the minimum accelerator-to-decelerator grid separation ratio, as 
in Fig. 21, but for a rather large accelerator-to-decelerator grid sepa- 
ration ratio of about 0.7. There is no obvious reason why this should 
occur. Figure 30 does show that accelerator-to-decelerator grid separa- 
tion ratio changes have only a minor effect on the minimum beam divergence 
angle magnitude. From Figs. 21 through 24 one would expect the ratio 
Jid/^s smaller effect as the net-to-total accelerating 

voltage ratio is increased. 

It is possible that at an accelerator-to-decelerator grid separation 
ratio of about 0.7 the equipotential surfaces between the screen and accel- 
erator grids adopt an optimum shape which results in the lowest off-axis 
velocity components for ejected ions and that this could effect the lowest 
divergence angles. However, it should be noted that an accelerator-to- 
decelerator grid separation ratio of 0.741 did not correspond to a minimum 

in decelerator grid impingement current. Decelerator grid impingement 

1 

current for the grid sets of Fig. 30 increased with accelerator-to- 
decelerator separation ratio in the same manner as those of Fig. 21. 

Further work is being done on the large screen-to-accelerator grid separa- 
tion ratio grid set of Fig. 30, to determine the minimum ion beam divergence 
angle possible with high net-to-total accelerating voltage ratios. 

Conclusion 

The use of three- grid accelerator systems in place of two- grid accel- 
erator systems offers significant improvements in ion beam focusing. These 
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improvements are greatest for low net- to- total accelerating voltage ratios 
and they should be negligible at a net- to- total accelerating voltage ratio 
of about 1.0; For a small screen- to-accelerator grid separation ratio, 
ion beam divergence decreases with decreasing accelerator- to-decelerator 
grid separation ratios and increasing net-to-total accelerating voltage 
ratios. This effect becomes small and is masked by second order focusing 
effects for a large screen- to-accelerator grid separation ratio. Decel era- 
tor grid impingement current increases sharply as either the accelerator- 
to-decelerator grid separation ratio is increased or the decel era tor hole 
diameter is reduced to the size of the accelerator hole diameter. Net-to- 
total accelerating voltage ratio variations do not affect decelerator grid 
impingement current as much as either of the parameters mentioned above. 
Initial results indicate that low beam divergence and low decelerator grid 
impingement currents for three-grid accelerator systems require that the 
decelerator grid be as close to the accelerator grid as possible and the 
decelerator hole diameter be about twenty percent larger than the accelerator 
hole diameter. Ion source beam current regulation is desirable when using 
three-grid accelerator systems to insure continued accelerator system opera- 
tion over the narrow beam current range at which the minimum beam divergence 
angle occurs. 



KEEPER ELIMINATION STUDY 


The major motivation for eliminating the keeper is the simplification 
associated with the removal of this electrode, i_ts_.power-supply and con- 
trols from“a thruster system. In order to determine whether its elimina- 
tion would degrade thruster performance, stability and/or operating and 
start up reliability a 15 cm SERT II thruster, equipped with high perveance 
dished grids, was tested in a mode where the keeper could be readily 
switched in or out electrically. This facilitated direct comparison of 
thruster start up and operation with and without the keeper. 

Test Description 

1 8 

A SERT II thruster equipped with high perveance small hole accelera- 
tor grids (67% open area screen and 24% open area accel ) was used for this 
test. It was modified so the cathode was isolated electrically from the 
thruster body. Further, it was connected to the required power supplies in 
the manner shown in Figure 31, through the switch "S", which was installed 
so the keeper could be isolated from or connected to its power supply, while 
the thruster was operating. High voltage tickler electrodes were not used 
to facilitate start up during any of the tests. During normal thruster 
operation the switch "S" in this figure is closed and the resistor "R" is 
an open circuit. Start up is achieved in this "conventional keeper" con- 
figuration by biasing the keeper a few hundred volts positive of the cathode 
to establish the keeper-cathode discharge. 

To test the effect of keeper elimination the switch "S" was opened so 
the keeper floated and, while it was still present in the thruster, it did 
not participate in the discharge process. Further the resistance of "R" 




FIGURE 31 
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was set at TOO so start up could be achieved in the manner suggested 

12 

by Kaufman, et al. With the thruster arranged in this configuration 
and the discharge extinguished, the entire thruster body, the screen grid 
and the cathode pole piece are biased to anode potential through the re- 
sistor "R". Anode potentials of 60 V were generally sufficient to effect 
start up although 100 V was required on one occasion. After start up the 
cathode couples to the thruster body, screen grid and cathode pole piece 
through the plasma thereby inducing a voltage drop across "R", The cur- 
rent to these thruster surfaces is sufficiently small that this voltage drop 
(the potential difference between the surfaces and the cathode) is of the 
order of several volts. The value of the resistor "R" was varied during 
the test and, while it would.be preferred to use a higher resistance to 
minimize losses in the resistor, it was found that with higher resistances 
the discharge would extinguish and not restart until the high voltage was 
turned off. With a 100 resistor installed the thruster ran stably and 
started reliably. 

Because only electrical alterations were required to shift between the 
two operating configurations, this change could be accomplished quickly and 
variations in start up and operating characteristics of the thruster due to 
propellant flow rate changes were eliminated. 

The Effect of the Keeper on Performance 

Figure 32 shows a typical comparison of the potentials of various 
thruster surfaces as a function of discharge current for the thruster 
operating 1) in the conventional way (at 0.3A keeper current) and 2) with 
the keeper disconnected. With the keeper disconnected it floats and 
provides an indication of the plasma potential in the cathode discharge 
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region. Figure 32 shows therefore that in the "keeper floating" configura- 
tion the keeper-to-anode potential difference is greater than it is in the 
"conventional keeper" configuration. At a given anode potential therefore 
primary electrons will acquire more energy as they pass into the main dis- 
charge region in the "keeper floating" configuration. In addition Figure 32 
shows that the thruster body operates positive of the cathode in the "keeper 
floating" configuration. One would expect a reduction in the ion sputtering 
damage of the screen grid, baffle and other thruster body-potential surfaces 
to accompany this change for the following reason. An increase in the po- 
tential of the thruster body relative to the anode potential implies a re- 
duction in the plasma-to-body potential difference. This in turn implies 
a reduction in the energy of the ions which strike the body- potential sur- 
faces and hence a reduction in the associated damage. 

Figures 33 and 34 show the effect of eliminating the keeper on the 
discharge loss-propellant utilization plots for the two total flow rate 
cases of 730 mA {full flow) and 435 mA-(throttled flow). The solid symbols 
on each curve correspond to operation at a 35 V discharge voltage. The 
performance penalty associated with the "floating keeper" configuration is 
observed to be negligibly small for both flow rates. It is interesting to 
note that part of the power going to the thruster bypasses the discharge 
chamber through the resistor "R" in the "floating keeper" configuration. 

This corresponded to discharge loss of 20 eV/ion at the solid symbol opera- 
ting point of Figure 33. If a switching circuit were employed which con- 
nected the thruster body to the anode when the discharge extinguished, in 
place of the resistor "R", this loss could presumably be eliminated. This 
should result in a lowering of the "floating keeper" curve by approximately 
20 eV/ion. 
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Doubly- to-singly charged ion density ratios were also measured using 
an E X B momentum analyzer during these tests. They suggested slightly 
higher doubly charged ion densities for the "floating keeper" configuration. 
Unfortunately the momentum analyzer was found to be giving irreproducible 
results shortly after these tests and the results are therefore considered 
unreliable and they are not presented here. 

Concl us ion 

The application of a 100 Si resistor between the anode and thruster 
body of a 15 cm mercury ion thruster provides sufficient coupling between 
these two components to facilitate thruster start up at normal anode volt- 
ages without the aid of a keeper. Operation in this configuration causes 
no significant degradation in thruster performance or stability and it 
results in a biasing of the thruster body and screen grid above cathode 
potential. This biasing should have the beneficial effect of reducing 
ion sputtering damage. 



MERCURY HOLLOW CATHODE STUDIES 
Dan Siegfried 

The operating temperature in and voltage drop across the orifice 

region of a hollow cathode are very important parameters which have been 

shown to affect the operational lifetime of the cathode. Prediction of 

this temperature and voltage drop as a function of cathode geometry and 

operating conditions requires a knowledge of the operating mechanism of 

the cathode and a model of energy transfer in the orifice region. At 

present the degree of understanding of these phenomena is limited. There 

have been numerous parametric studieb^^’ of hollow cathodes and 
22 

Bessling has proposed a theoretical model of the cathode, but verifica- 
tion of the model as well as an understanding of the basic physical proces- 

23 

ses taking place in the cathode is still incomplete. Finally, Krishna%\ 
has tested and modelled the physical processes occurring in a large 
diameter argon hollow cathode driven by a pulse line. The applicability 
of his collision model to this smaller mercury device remains to be verified. 
The objective of this study has been therefore to gain a better insight 
into the physical processes of hollow cathode operation and the applicability 
of these models by using Langmuir probes to map plasma properties upstream 
and downstream of the cathode orifice plate. A comparison of these experi- 
mental results with those predicted by the models in describing the cathode 
phenomena observed has also been initiated. 

Apparatus and Procedure 

All tests were conducted on a hollow cathode composed of a 2% thoriated 
tungsten orifice plate electron-beam welded to a 6.35 mm diameter tantalum 
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tube. The orifice plate was 1.22 mtn thick with a 0.75 mm diameter orifice 
on its centerline. The cathode contained a rolled tantalum-foil insert 
coated with chemical R-500 and was mounted in the operating structure shown 
in the cutaway sketch of Figure 35. Langmuir probes were mounted upstream 
and downstream of the orifice on the centerline of the assembly. They could 
be moved axially over a 0 to 25 mm distance from the orifice plate during 
cathode operation. The upstream probe entered the rear of the cathode 
fixture through a feedthrough fixed to the cathode body. In order to mini- 
mize leakage of mercury vapor past the probe the feedthrough was made 
approximately 8 probe diameters long and the probe- to-feedthrough clearance 
was held to about .05 mm. 

Probes were constructed using quartz tubing drawn to a close fit 
around the tungsten electrode. Based on the initial testing, it was de- 
termined that the upstream probe electrode should be .25 mm diameter wire 
flush with the end of the quartz insulator and the downstream electrode 
should be .63 mm wire exposed for a length of approximately 1.0 mm. Using 
these probes, voltage-current traces were recorded at various positions up- 
stream and downstream of the orifice plate. These measurements were recorded 
for anode currents ranging from .5 to 9.0 amps and mercury flow rates of 
100 and 150 mA. Both spot and plume modes of operation were observed at 
these flow rates. 

After a portion of the testing had been completed it was found that 
improved resolution could be obtained using much finer electrodes and 
drawing the quartz tubing to ''^.5 mm diameter for 5 mm behind the exposed 
portion of the electrode. For these tests the exposed electrode wire was 
.075 mm diameter by .25 mm long for the upstream probe and .125 mm diameter 
by 1.0 mm long for the downstream probe. Using these probes, data were 



HOLLOW CATHODE TEST CONFIGURATION 





-72- 


recorded at a 100 mA flow rate for both 2.0 A (plume mode) and 6.0 A 

(spot mode) discharge currents. For all tests the keeper current was held 

constant at 0.3 A. Langmuir probe traces were analyzed in all cases using 

24 

the technique proposed by Beattie. 

Experimental Results 

Discharge current- voltage plots for the two flow rates investigated 
are shown in Figure 36. The regions of spot and plume mode operation are 
also indicated. The plasma properties pertaining to the currents and 
operating modes of Figure 36 are shown in Figures 37 through 40.- These data 
were all obtained from measurements with the larger diameter probes men- 
tioned under "Apparatus and Procedure." The curves are notable in that they 
show magnitudes of plasma properties, their variation with position and 
changes in the plasma associated with transitions between the spot and plume 
modes. For example it can be seen that in the spot mode at both flow rates 
the electron density profile is independent of the discharge current. In 
the plume mode however reductions in discharge current induce a corresponding 
reduction in electron density. The electron temperature and most plasma 
potential profiles also appear to be more strongly affected by discharge 
current in the plume mode. Downstream of the orifice in the highly luminious 
plume mode, both the plasma potential and the electron temperature increase 
with increasing discharge current. When the discharge switches from the 
plume to the spot mode, the curves for these properties drop significantly 
and change shape. This drop is illustrated by comparing the 3.0 A and the 
7.0 A curves for the 100 mA flow rate case (Fig. 37) or the 5.0 A and 2.0 A 
curves for the 150 mA flow (Fig. 39). The 7.0 A and 5.0 A currents repre- 
sented the lowest current levels for stable operation in the spot mode (see 
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FIGURE 37 






UPSTREAM 


DOWNSTREAM 


DISTANCE (mm) 


CATHODE PLASMA POTENTIAL AND ELECTRON TEMPERATURE PROFILES 

150 mA FLOW RATE 

FIGURE 39 
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Figure 36) and visually corresponded to a very intense spot at the orifice 
with very little luminosity present away from this spot. As the discharge 
current is increased above this value, the plasma potential and electron 
temperature (Figs. 37 and 39) again increase in what appeared to be direct 
proportion to the intensity of the radiation from the spot. 

It should be noted that it was not possible to probe the region from 
about 4 mm upstream of the orifice to 1 or 2 mm downstream for either mode 
due to the very high electron density in that region, so the curves can 
only be estimated there. The data scatter on the electron temperature is 
about .25 eV and on the plasma potential about one volt except for the 
0.5 A and 1.0 A cases where it is higher. The data scatter for the electron 
density is less than an order of magnitude on the semi -log, scale plots. 

The fact that the plasma potential increases from a location some 15 mm 
upstream of the orifice and for some distance downstream of the orifice, as 
shown in Figures 37 and 39, is significant. The electron collision rate 
must be very high throughout the region to sustain this potential gradient. 
Because electron temperature is not increasing proportionately over the same 
region it appears that energy transferred to the electrons as they pass 
through this potential gradient must be lost through the collisions or by 
radiation. 

Plasma property results measured using the smaller diameter Langmuir 
probes mentioned under "Apparatus and Procedure" are shown as solid lines 
and data points in Figures 41 and 42 for the plume mode and Figures 43 and 
44 for the spot mode. These plots are very similar to the corresponding 
ones made with the larger probes. Figures 37 through 40, except that the 
steep gradients in the plasma potential and electron temperature in the 
region 5 to 10 mm downstream of the orifice in the plume mode were not 
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FIGURE 43 
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readily discernable from the large probe data. The presence of this spike 
was, however, suggested in the earlier data by a few stray data points in 
this region that were originally attributed to data scatter. The better 
resolution of the smaller probe allowed the more accurate mapping of this 
region. 

The steep electron temperature and plasma potential gradients in 
Figure 41 apparently correspond to an electron free-fall region where the 
electrons are free to be accelerated through a potential difference of about 
20 V without significant collisional interaction. That this is occurring 
is born out both by visual observation of the discharge and analysis of the 
probe traces. Observation in the plume mode showed a dark band followed by 
a region of more intense luminosity as shown in Figure 41. At a discharge 
current of 2.0 amps this band was located about 5 mm downstream of the 
orifice. Neither the dark band nor the spikes in plasma potential and 
electron temperature were present in the spot mode. It should be noted 
that the intense spot was present in both modes (Figures 41 and 43). It 
was slightly larger and noticeably more intense in the spot mode. 

Computer analysis of the Langmuir probe traces indicated that near the 
peak in the plasma potential in the plume mode curve (Fig. 41), about 30% 
of the electrons were monoenergetic (primary) and that the energy of this 
electron group was about 19 eV. This agrees well with the ''^20 V rise in 
plasma potential which takes place between 1 and 6 mm downstream of the 
orifice. Electrons passing through the free-fall region would have suffic- 
ient energy to cause the luminosity seen immediately downstream of the dark 
"free-fall" region where inelastic exitation reactions are presumably 
occurring. The presence of collisions is indicated in the plasma property 
data by the decrease in the primary fraction and the eventual randomization 
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of electron energies toward a Maxwellian distribution of about 2 eV as 
the electrons diffuse away from the free-fall region (dark band) shown in 
Figure 41 . 

The solid symbols on the electron temperature and plasma potential 
curves are questionable data points (Figure 41). The computer analysis 
of the associated Langmuir probe traces suggested no primary electrons 
were present, and then indicated rather large errors were brought on by 
a pure Maxwellian curve fit. This suggests the Maxwell ian-plus-mono- 
energetic distribution does not fit well at these points. They are, however, 
considered approximately correct and they do fit with the other data. It 
should also be noted here that the banded appearance indicated in Figure 41 
and discussed above was very faint and without sharply defined boundaries. 
This band has, however, since been observed under other operating conditions 
with .38 mm diameter orifice cathode where it was much more distinct. The 
dark band was observed to be very close to the luminous spot at low dis- 
charge current levels in the plume mode and to move 5 to 10 mm downstream 
as the current was increased. When the discharge changed to the spot mode 
this dark band seemed to disappear abruptly. 

Two physical processes important to hollow cathode design are electron 
emission and heat transfer rate to the orifice. If it were possible to pre- 
dict the plasma conditions at the orifice, then the calculation of both of 
these would be facilitated. Two approaches to the heat transfer rate calcu- 
lation might be taken. In one continuum flow through the orifice would be 

assumed and a convective heat transfer coefficient based on local plasma 

22 

properties would be employed. Bessling suggests another model using a 
particle approach. Either approach requires assumptions which have not 
previously been amenable to verification. In the continuum approach an 
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important assumption is thermal equilibrium and the applicability of the 
Saha equation. Bessling's model is highly dependent on- the value of an 
assumed ionization fraction in the orifice region— a quantity which is 
also important in Bessling's modeling of the electron emission process. 

The results of Langmuir probe analysis would be most useful therefore if 
they could be obtained in the region of the cathode of primary interest to 
these models, namely the orifice region. Unfortunately reliable probe data 
could not be obtained within the region from 5 mm upstream to 2 mm down- 
stream of the orifice without burning up the probe. Further, interpolated 
data are not reliable because of the large changes in density which appear 

to occur through this region. Through the application of an expanded Saha 
. 22 

equation and the designation of reasonable temperatures and pressures in 
the region upstream of the orifice, however, some insight into the required 
variation in plasma properties in this region can be obtained. 

In the development of his hollow cathode model Bessling cites the 
expanded Saha equation mentioned above in which the electrons and neutral 
atoms/ions though presummably in equilibrium with themselves are at two 
different temperatures and therefore not in equilibrium with each other. 
This equation is given as 
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where Tg is the neutral atom and ion temperature (eV), T^^ is the electron 
temperature (eV), a is the ionization fraction, P is the total (electron 
plus heavy particle) pressure, m^ is the electron mass, h is Planck's 
constant, e is the electron charge, is the ratio of internal partition 

functions for the singly ionized and neutral ground states (2 for mercury), 
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i{i. is the- ionization potential, and G is the ratio Tp/T . The total 
1 rnx 

pressure is given by 

P - "e "g 

where n^ and are the electron and heavy particle (atom plus ion) 
densities respectively. The Langmuir probe data provide values for T 

lltA 

as a function of axial location. If and P were known then the ioniza- 
tion fraction a and electron density n could be calculated from 
Equations 6 and 7 above at each axial location. The pressure P can be 
estimated by assuming that the pressure in the cathode upstream of the 
orifice is uniform and is sufficient to insure choked flow through the 
orifice at the flow rate being maintained for the tests. Then, if energy 
addition in this region is neglected, the stagnation pressure is a function 
only of gas temperature and known physical parameters according to the 
equation , 



where m is the neutral propellant flow rate, D is the orifice diameter, 
k is the ratio of specific heats for the vapor (assumed to be near 5/3 
although the results are not sensitive to its value), and R is the gas 
constant. The temperature Tg is an unknown but bounds can be placed on 
it from consideration of the measured temperatures prevailing in the cathode. 
Its value should lie between the local electron temperature and some 

* The partial pressure of the electrons increases while that of the neutrals 
and ions decreases as one proceeds downstream toward the orifice. This 
is assumed to occur in such a way that the sum of these two remains con- 
stant. In reality the total pressure is probably fairly uniform until 
the immediate vicinity of the orifice is approached and a slight decrease 
in pressure is observed. 
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temperature characteristic of the internal cathode surfaces. The maximum 
electron temperature obtained in the cathode was about 1.0 eV and cathode 
surface temperatures above 150°C (a- 0.036 eV) have been found necessary to 
prevent mercury condensation. Using these two values as bounds the pres- 

25 

sure must lie between .75 and 4.0 torr. Hollow cathode temperature studies 
have shown that the average temperature near the orifice is closer to 

.1 eV (880°C) so this was' selected as the most probable temperature for 

the internal surfaces. Using this temperature and the internal pressure of 
1.25 torr calculated on the basis of it, the electron density {n ) and 
ionization fraction (a) profiles were calculated for the values of Maxwellian 
electron temperature shown in Figure 43 as a typical case of spot mode 

operation. The corresponding values of electron density are shown in 

Figure 44 plotted as broken lines along with the experimental data (solid 
lines). Values of a are plotted in Figure 45. To determine an upper bound 
on the density and ionization fraction they were also calculated for the 
case of a 1.25 torr pressure and a gas temperature equal to the local elec- 
tron temperature. These are also plotted on Figures 44 and 45 as center 
lines. 

It was apparent that similar results would not be obtained for the 
electron temperature profile for the plume mode shown in Figure 41 as the 
solid line because the line was too flat. In this case the temperature 
profile shown by the dotted line was assumed and the same calculations dis- 
cussed above resulted in the theoretical curves indicated on Figures 42 and 
46. Although the assumed (dotted) temperature profile deviates from the 
best fit for the data points for upstream of the orifice it is seen to agree 
fairly well with the data points within 10 mm of it. An error bar shown at 
the 7.5 mm position on the calculated curves for constant gas temperature 
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(dashed lines) in Figures 42, 44, 45 and 46 indicates the range of density 
variations calculated for the cases where the neutral /ion temperature (T„) 
assumes the maximum probable range of internal cathode surface temperatures. 

The results of the above calculations are somewhat sensitive to the 
simple assumptions made regarding the magnitude and shape of the gas tem- 
perature profile and the total pressure profile in the cathode chamber. 
However, for such a simple model, the agreement with the experimental data 
is encouraging, and allowing Tg to vary should facilitate even better agree- 
ment. For example the close proximaty of the solid and dashed lines at 
locations more than 5 mm upstream of the orifice, suggests' the ion/neutral 
temperature is probably about. equal to the internal cathode surface tem- 
perature in this region. In this region then the ions and neutrals are not 
in equilibrium with the electrons. In the region zero to 5 mm upstream of 
the orifice there is no straightforward way to tell if the neutrals and 
ions are in equilibrium with the electrons. At the orifice itself however 
in the cases of the lower temperature Tg (dashed lines on Figs. 45 and 46) 
Equation (6) requires that the ionization fraction come infinitesmally close 
to unity. In the case where the temperature Tg is equal to the electron 
temperature (center lines on Figs. 45 and 46) both the shape of the curves 
at the approach to the orifice and the magnitudes of the ionization fractions 
are more reasonable (80 to 90%). This suggests that the atom/ion tempera- 
ture probably increases toward the electron temperature as the orifice is 
approached. 

The results presented in Figures 41 to 46 can now be used to gain some 
insight into emission mechanisms in the' hollow cathode. Bessling concluded 
on the basis of his model that the dominant emission mechanism was thermionic- 
field emission in the ori-fice region. In order to evaluate this, note that 
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in a plasma where' the electrons and ions are at two different temperatures 
the Debye length, Xp, is given by 



where is the permittivity and n^ is the ion density apd is equal to the 
electron density. Figures 42 and 44 suggest that at most the electron 
density is 1.2 x 10^^ cm"^. Even assuming that the gas temperature and 
electron temperature differ considerably say Tp = .1 eV and T = .8 eV 
the Debye length would be no smaller than 6.33 x 10“® cm. The plasma 
potential in the orifice region is about 10 v. These two results combine 
to suggest an electric field strength at the orifice wall of 

E = f- = 1.58 X 106 v/cm. (10) 

^D 

This is an order of magnitude below the_field strength of 10^ V/cm 

22 

suggested by Bessling as necessary for the electric field to contribute 
significantly to the emission. Without accounting for the possibility of 
locally high fields due to surface irregularities, it appears that 
thermionic field emission is insignificant in the orifice region. 

Cathode-Insert Interaction Studies 

One additional test was conducted on a hollow cathode like the one 
shown in Figure 35 which warrants some discussion. In this test the 
rolled foil insert shown in the Figure was replaced by a 1.9 mm dia. 
cylindrical porous tungsten insert which had been impregnated with low 
work function material. (Semi con type 84-S). In addition the Langmuir 
probes were not present and the insert was attached to a wire which passed 
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through the feedthrough to facilitate axial movement of the insert during 
cathode operation. A new cathode, which had never been exposed to low 
work function oxides and which had a 0.76 mm dia orifice was selected for 
the test. The upstream plate shown in Figure 35 through which the insert 
moving wire was fed was made of quartz for this test. The quartz window 
facilitated direct observation of the discharge as well as temperature 
measurement of the interior surface of the cathode orifice plate using 
an optical micro- pyrometer. 

After the mercury flow rate had been stabilized' into the preheated 
cathode, the discharge was started with the downstream end of the insert 
positioned 1.75 cm upstream of the orifice plate. The cathode was 
operated for a period of several .minutes at a keeper current of 0.3 A. and 
a discharge current of a few amperes to allow the discharge to stabilize. 

The discharge current was then increased to 5 A. and the porous tungsten 
insert was moved toward the cathode orifice plate with the trajectory 
shown in the upper plot of Figure 47. It was^then pulled back to the 
1.75 cm position where it was held for two minutes before it was inserted 
and retracted rapidly as shown in Figure 47. The time variation in anode 
voltage and orifice plate temperature (measured at the edge of the orifice 
itself) which accompanied this insert movement are also shown in the 
figure. Although some of the times associated with the figure were 
estimated, the time scale is considered accurate to within t 15 sec. and 
data shown as being collected at the same time were collected with a few 
seconds of each other. 

From Figure 47 it can be seen that at the initial operating condition 
(a new cathode and the insert 1.75 cm from the orifice plate) the orifice 
plate temperature and voltage were both high 1400°C and 15 V respectively). 
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Mo Vi ng the insert downstream is seen to reduce both of these quantities. 
Moving the, insert back upstream is observed to effect a gradual increase 
in temperature and voltage apparently as a result of gradual depletion of 
the low work function material in the vicinity of the orifice. The corre- 
lation between anode voltage and temperature is striking. The temperatures 
shown in Figure 47 were measured at the edge of the orifice. Temperatures 
at the intersection between the orifice plate and cathode tube were observed 
to be about 100°C lower than these plotted values. After this test, the 
cathode was replaced with one which again had not come into contact with 
any low work function material. This cathode also had a 0.76 mm dia orifice 
but it had a dummy insert (no low work function material). It was observed 
to operate at orifice plate temperatures of about 1600°C and discharg volt- 
ages of 13 volts when the discharge current was 5 A. Contrary to the case 
with the insert-equipped cathode^ where temperature gradients were observed 
on the orifice plate, the entire orifice plate of this second cathode appeared 
to be at a uniform temperature. Further, the luminosity from excited mercury 
atoms could be seen throughout the entire cathode support structure and into 
the 1/8” propellant feed tube (Figure 35). When a low work function insert 
was present, luminosity was not observed up stream of the insert. 

Concl us ions 

Plasma potentials and electron temperatures appear to increase through 
the orifice of a cathode operating in either the spot or plume modes. 

Electron density profiles are very peaked in the vicinity of the orifice. 

The electron densities at a given location increase with discharge current 
in the plume mode and tend to remain constant with increases in discharge 
current in the spot mode. 
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Ups tream of the orifice th'e measured electron temperatures and 
densities correlate well with calculations based on the expanded Saha 
equation if the neutral/ion temperature is assumed to be close to the 
wall temperature and the pressure in the cavity is assumed to be constant. 
The ori.fice itself appears on the other .hand to be a region of a high 
electron collision frequency where thermal equilibrium between electrons, 
ions and atoms is approached. The results suggest that the ionization 
fraction in the orifice is high. They do not, however, suggest that 
thermionic-field emission is the dominant emission mechanism for the 
hollow cathode. 

The movement of a cathode insert into a previously virgin cathode 
suggests that the nature of the discharge changes markedly as the avail- 
ability of low work function materials in the vicinity of the orifice 
changes. When this material is not available discharge voltages and 
cathode temperatures increase greatly, and electrons tend to be drawn 
from surfaces far from the cathode orifice in an apparent effort to 
satisfy the requirement for the emission current demanded by the dis- 
charge. 



8 cm THRUSTER TESTS 


An 8 cm mercury ion thruster was assembled from component parts of 

the SIT 8 thruster^®’ and tested in the large {1.2 m dia. x 4.6 m long) 

vacuum facility at Colorado State University. The principal purpose of 

these tests was to determine if the doubly charged ion production model 
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proposed by Peters could be applied to discharge chambers as small as 
this 8 cm dia. one. 

Apparatus and Test Procedure 

The thruster was operated using two different sets of grids both of 
which employed 67% open area screens. The SHAG optics set had a 19% open 
area accel grid while the other one had a 44% open area accel grid. 
Verification of the doubly-charged ion production model required simul- 
taneous measurement of plasma properties within the thruster and the doubly 

and singly-charged ion density profiles in the thrust beam. The required 

29 

plasma properties were determined from Langmuir probe data collected at 
a total of at least 16 points within the discharge chamber at each opera- 
ting condition, and each of these data sets analyzed .using a computer 

24 

program developed for this purpose. The doubly and singly-charged ion 

current density profiles were measured in the ion beam of the thruster 

30 

using an E x B momentum analyzer. 

The thruster was operated at various flow rate and discharge voltage/ 
current conditions using the two grid sets and the aforementioned data 
were collected. Plasma property data were next volume averaged over the 
primary electron region within which the bulk of the ionization is 
presumed to occur. ' This volume averaging process, which is described in 
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detail by Peters, also enables one to determine the ratios of singly 
and doubly-charged ion fluxes within the primary electron volume to those 
on its surface. These ratios are necessary for the calculation of doubly 
and singly-charged ion loss rates in the model. 

Results 

During the test sequence operation at the nominal, condition was 
attempted to check if the source was performing as similar units had in 
other laboratories. Using the 19^ open area accelerator grid (SHAG) it 
ran stably and the operating variables shown in Table I were measured. 

The nominal values operating variable as supplied by NASA for the SIT 8 
unit are also listed for comparison. 


Table I 

8 cm Thruster Nominal Operating Condition 


Variable 

Nominal Value 

As Measured 

Anode- to-Neut. Potential 

1205 V 

1200 V 

Accelerator Voltage 

-300 V 

-300 V 

Accelerator Drain Current 

0.17 mA 

0.4 mA 

Discharge Neutral Flow Rate 

91.4 mA 

90 mA 

Beam Current 

72 mA 

73 mA 

Discharge Current 

0.455 A 

0.60 A 

Discharge Voltage 

35 V 

38.9 V 

Keeper Current 

0.1 A 

0.1 A 

Keeper Voltage 

6.0 V 

10.0 V 

Cathode Heater Current 

0 

0 

Neut. Keeper Current 

0.5 A 

0.5 A 

Neut. Keeper Voltage 

13.0 V 

20 V 


The most striking differences between the nominal and measured values in 
this table are related to the higher keeper voltage measured on both the 
main and neutralizer keepers. The high main keeper voltage necessitated 
in turn an increase in the discharge voltage so the discharge-keeper 
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voltage difference could be maintained at a specified value of about 29 v. 

The reason for the higher keeper voltages is not certain^but it is believed 
to be related to possible contamination of the main and neutralizer cathodes 
an/or depletion of the low work function material which should be present 
in them. This is substantiated by the fact that in the first tests of 
these cathodes^ main and neutralizer keeper voltages of 11 and 21 v res- 
pectively were observed. A subsequent attempt to add chemical R-500 
through the cathode orifices reduced these voltages to about 8 and 17 v 
respectively. After a few hours of operation however these voltages rose 
to the "as measured" values given in Table I. 

The "as measured" discharge current in Table I is also observed to 
be high. While this can be explained in part by the higher "as measured" 
beam current and lower discharge neutral flow rate, the discharge current 
is still considered high. These discrepancies were however considered 
sufficiently small so that the main objective of the tests (verification 
of the doubly-charged ion production/loss model ) could be pursued. 

The operating conditions at which 8 cm thruster data were collected 
for the verification of the doubly-charged ion model are given in the 
first five rows of Table II. The next row of this table indicates the 
volume-to-surface area ratio for the primary electron region of the 8 cm 
thruster as determined from iron filings maps is 0.92 cm. Following this 
quantity, the average plasma properties as determined from Langmuir probe 
traces; which had been analyzed, appropriately weighted and volume averaged; 
are listed. The uniformity factors which are required in the analysis are 
also listed. They are determined from measured properties and represent 
the primary electron region volume-to-surface area flux ratio for singly 
(F^) and doubly charged ions. Following this, the ratio of doubly-to- 
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Table II 


8 cm Thruster Calculated Densities and 
Production Routes for Various Operating Conditions 


Screen/Accel Grid Open Area Percentages 

Discharge Current (A) 

Discharge Voltage (V) 

Beam Current (inA) 

Mass Flow Rate (mA) 

67/44 

.71 

42.7 

66 

86 

67/44 

.26 

39.2 

44 

90 

67/19 

.63 

,38.6 

73 

90 

67/19 

.86 

41.0 

77 

91 

67/19 

.66 

40.0 

72 

88 . 

67/19 

.60 
39 
73 
90 • 

67/19 

_. .70 

39 
;■ 75 
91 

67/19 

.38 

• 37,6 
56 
. 91 


^Plasma Volume- to-Surface Area Ratio (cm) 

.92 

.92 

.92 

.92 

.92 

.92 

.92 

.92 


Avg. Maxwellian Electron Temperature (eV) 

13.0 

5.2 

7.6 

. 12.6 

9.7 

7.7 

8.8 

4.2 


Avg. Primary-to-Haxwellian Elec, Density Ratio 

0.19 

0.052 

0.11 

0.087 

0.089 

0.078 

0.061 

0.029 


Avg. Primary Electron Energy (eV) 

29.0 

24.7 

27.3 

23.2 

19.5 

21.3 

19.4 

20.0 


Avg. Electron Density (cm"3 x 10"'^®) 

,9-2 

8.5 

12.6 

12.5 

11.0 

11.8 

12.9 

■ 14.4 


Uniformity Factor F+ 

2.11 

1.77 

1.93 

i;8 

1.85 

1.81 

1.88 

1.53 


Uniformity Factor, F^^ 

1.99 

1.99 

1.92 

2.02 

2.67 

2.63 

2.74 

1.56 

ip 

CD 

1 

Measured Doubly- to-Singly Charged Ion 
Current Percentage 

7.96 

2.0 

5.4 

15.4 

9.5 

4.6 

6.2 

1.5 

lieutral Ground State Atoms 

' .70 

.64 

.62 

.64 

.63 

.61 

.61 

.55 


Neutral Metastable State Atoms 

.033 

.084 

.065 

.039 

.055 

.068 

.063 

.11 


Neutral Resonance State Atoms 

.21 

.26 

.26 

.26 

.26 

.28 

.27 

' .32 


Singly Charged Ground State Ions 

.060 

.021 

,055 

.067 

.049 

.044 

.054 

.021 


Singly Charged Metastable Ions 

.001 

.000 

.001 

.001 

.002 

.002 

.002 

.001 


Doubly Charged Ground State Ions 

.002 

.000 

.001 

,002 

.001 

.001 

.001 

.001 

O o 

. Neutral Ground State 

.63(.27)* 

.47(.32) 

.50(.33) 

.56(.12) 

.52(,13) 

.48(.18) 

.49(.ll) 

.34(.24) 

Neutral Metastabie States 

.066(.25) 

.19(.21) ■ 

.13(.26) 

.077(.ll) 

.llK.12) 

.14(.15) 

.13(.09) 

.27(.14) 

Es 

Neutral Resonance States 

.31(.27) 

.34(.28) 

.37(.31) 

.37(.ll) 

.36(.12) 

.38(.17) 

.38(.09) 

.39(.18) 

[^Neutral Ground State 

.25(.037) 

.093(0.) 


.23(0.) 

.22(0.) 

.14(0.) 

.17(0.) 

.053(0.) 


Neutral Metastable State 

.012(.037) 

.012(0.) 

.010(.017) 

.014(0.) 

.019(0.) 

.016(0.) 

.018(0.) 

.011(0.) 

c? 

Neutral Resonance State 

.075(.037) 

.037(0.) 


.094(0.) 

.089(0.) 

.065(0.) 

.078(0.) 

.030(0.) 

^ Q 

Singly Charged Ionic State 


.85(.48) 

.84(.41) 

.65(.13) 

.67(.092) 

•77(.24) 

.72(.07) 

.90(.43) 

^Singly Charged Metastable State 

;.005(.22) 



.008(.026) 

.009(.034) 

.009(.03) 

.012(.003) 

1 .003(.05) 


* Numbers in parenthesis are the fractions of the indicated interactions effected by primary electrons. 
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singly-charged ion current in the beanijas determined from mass spectrometer, 
data are listed. Dividing each of the numbers in this row by 2 yields 
the average doubly- to-singly charged ion density ratio in the discharge 
chamber expressed as a percentage. 

The section in Table II titled “Calculated Fractions of lons/Atoms 

in Various States" lists the normalized densities of the states included 

in the analysis and predicted by the model. The normalized density of a 

given specie is defined here as the specie density divided by the total 

heavy particle density. The sum of the normalized densities for any 

thruster operating condition should therefore equal unity. Table II shows, 

for example, that the 8 cm diameter thruster operating with the 67/44 grid 

set at a 42.7 v discharge voltage would be predicted to have 70% neutral 

ground state atoms, 21% neutral resonance state atoms, 6.0% singly-charged 

ground state ions and .2% doubly-charged ground state ions. The data in 

this section of the table are observed to be quite similar to those re- 
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ported by Peters for the 15 cm and 30 cm thrusters. 

The last section in Table II shows the calculated production rates 
for singly and doubly-charged ions through the various intermediate states 
listed. These production rates have been normalized by the total produc- 
tion rate of the specie indicated. The fraction of the associated inter- 
actions effected by the primary electrons is indicated in parenthesis. 

For example, at the thruster's 42.7 V, 67/44 grid operating point, 63% of 
the single ions are produced as a result of electron interaction with 
neutral ground state atoms and 31% are a result of electron bombardment of 
neutral resonance state atoms. The neutral ground state-to-single ionic 
state interactions were induced by primary electrons 27% of the time and 
by Maxwellian electrons the remainder (73%) of the time. Again these 
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values are relatively similar to those reported by Peters for the 15 cm 
diameter thruster. Comparison of the data of Table IT with similar data 
for the 15 cm and 30 cm thruster reveals the smaller the thruster, the 
larger the electron temperature and the greater the fraction of the inter- 
actions that are induced by Maxwellian electrons. In all three thrusters 
doubly-charged ions are produced chiefly through the intermediate singly- 
ionized ground state. 

The ratio of doubly-to-singly charged ion densities in the discharge 
chamber is the most important quantity generated by the analysis which pro- 
duced the data of Table II, The extent to which this ratio agrees with the 
measured ratio is'shown in Figure 48 for each of the operating points cor- 
responding a column of Table II. Most of the data on this figure are ob- 
served to fall close to the 45° line corresponding to perfect correlation. 

One point corresponding to operation at the highest doubly charged ion 
density ratio is observed to lie significantly outside of the t 20% error 
band, but in general the correlation appears to be good. 

In discussing the errors in this work it is necessary to point out 
that the uncertainty in plasma properties, obtained from analysis of Langmuir 
probe traces, was observed to be significantly greater in the 8 cm thruster 
than it had been in the 15 cm one. It is believed that this increased un- 
certainty occurs because the 8 cm thruster plasma departs substantially 
from the Maxwell ian-plus-monoenergetic distribution function that is an 
integral part of the probe trace analysis procedure. 

Electron Acceleration Through the Baffle Aperture 

A Langmuir probe was installed within the 8 cm thruster so it could be 
moved axially through the baffle aperture during thruster operation in the 
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manner suggested by the sketch at the bottom of Figure 49. Measurements 
were made using this probe with the thruster operating at the conditions 
indicated on Figure 49 and with a open area accel grid. The variation 
in the properties of the cathode and main discharge plasmas through the 
baffle aperture are also shown in Figure 49 It is considered significant 
that electron temperature and "plasma potential show a rather steep increase 
10 to 12 mm upstream of the aperture and primary electrons are observed 
some 6 mm upstream of the aperture. In fact it appears from the primary 
electron energy and plasma potential curves that the acceleration process 
has been -completed by the time the electrons reach a position 6 mm upstream 
of the baffle. 

These high^^energy electrons located upstream of the baffle can enter 
] 

into reactions producing excited atoms, and singly and doubly-charged ions 
which will probably not escape from the cathode discharge region. These 
represent unnecessary energy losses which could be reduced if the electron 
acceleration process could be effected “in the more immediate vicinity of 
the baffle aperture. It is believed this could be accomplished if the 
magnetic field lines coming from the cathode pole piece could be concen- 
trated near its downstream edge. 

Conclusions 

The model for doubly-charged ion production and loss developed by 
Peters can be applied to the 8 cm thruster with reasonable accuracy 
{y i 20%). The bulk of the doubly charged ions in an 8 cm dia thruster 
are produced from singly charged ions ionized by electrons from the high 
energy tail of the Maxwellian group. 

Primary electron acceleration in the 8 cm thruster occurs significantly 
upstream of the baffle aperture. 
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SCREEN ANODE THRUSTER 


Glen Longhurst 
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Life tests conducted on the 30 cm mercury bombardment thruster 
have shown screen grid erosion becomes life limiting under some thruster 
operating conditions. This erosion is cause by ions which are accelerated 
through the plasma-to-screen grid potential difference and impact on the 
screen grid. The bulk of the damage appears to be caused by doubly-charged 
ions, which acquire twice the energy of singly charged ions, as they are . 
accelerated through a given sheath voltage drop. The sputtering yield S 
of a surface (atoms ejected per incident ion) in the energy range of in- 
terest here is given by 

S a (E-E^}3 ^ (11 ) ■ 

where E is the incident energy of the ion (typically of order 80 eV for a 
doubly-charged ion) and E^ is the sputtering threshold energy of the target 
material (approximately 16 eV for molybdenum). The cubic dependence of 
sputtering yield on energy suggests that even modest reductions in the in- 
cident ion energy will reduce the sputtering damage substantially. One 
method of reducing this ion energy would be to raise the potential of the 
screen grid to a value closer to plasma potential. Small increases in 
screen grid potential, such as those described previously in the "Keeper 
Elimination" section of this report, can be effected without changing the 
design chamber design and the thruster will still perform well. When the 
screen grid potential is increased to the point where primary electrons 
will have sufficient energy to' pass from the plasma to the screen grid, 
however, the grid becomes an anode. In this case the thruster must be re- 
designed so its magnetic field prevents direct primary electron collection 
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on the screen grid. 

One thruster design in which the screen grid is shielded by the 
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magnetic field is the radial field thruster design. A sketch of one 
configuration of this thruster with an iron filings map of its magnetic 
field shown as Figure 50 illustrates the salient features that make it 
suitable for this application. Magnetic field lines are observed to pass 
essentially radially from the center pole piece to the outer pole piece 
so primary electrons emitted through the holes in the center pole piece 
tend to be confined between the dotted critical field- lines. As the 
figure suggests therefore, primary electrons do not have direct access to 
either the screen grid or the upstream anode and both of these surfaces , 
can be biased positive without depleting the primary electron component 
of the discharge plasma. The thruster is designed so the radial and axial 
electromagnet currents can be varied independently and so the upstream 
anode can be moved axially during thruster operation to facilitate thruster 
performance optimization. The thruster employed for these tests was 15 cm 
in diameter and employed a 67% open area screen grid and 54% open area 
accelerator grid. 

Effect of Biasing Screen Grid Toward Anode Potential 

Numerous changes were made in the discharge chamber design in an 
effort to reduce the discharge losses and increase the propellant utiliza- 
tion associated with its operation. Design deficiencies which resulted 
in poor performance in early tests, were associated with 1) the establish- 
ment of a high impedance at the point of electron acceleration into the 
main discharge through the holes in the center pole piece and 2) the 
production of a primary electron region of insufficient axial depth and 
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uniformity to produce a high propellant utilization. While these de- 
ficiencies have not been corrected to the point where the performance is 
near that reported by Poeschel or observed in other 15 cm thruster 
designs, it is considered sufficiently good to permit evaluation of the 
concept of screen grid biasing. Figure 51 shows the best performance 

curve achieved with the radial field thruster compared to that observed 
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with the SERT II thruster operating with the same grid set and at the 
same propellant flow rate. The solid symbols on both curves represent 
operation at a 37 V discharge voltage. The effect of biasing the screen 
grid toward this anode voltage of 37 V while holding the upstream anode volt- 
age constant at 37 V and total discharge current constant at 4.1 A is illus- 
trated in Figure 52. Ion beam current (and propellant utilization), dis- 
charge losses and the fraction of the discharge current carried by the 
screen grid are shown plotted against the screen grid bias above cathode 
potential. The curves show discharge losses increase and propellant utili- 
zation decreases as the screen grid bias is increased and the fraction of 
the discharge current collected by the screen grid increases. While this 
is undesirable, the data do show that this decrease in propellant utiliza- 
tion and increase in discharge losses is quite small until the screen grid 
bias is about 30 V (7 V below anode potential). This data suggests that 
the screen grid could be biased 10-15 V above cathode .potential where 
sputter erosion damage would be greatly reduced without the occurrence 
of substantial performance degradations. Additional improvements in per- 
formance are possible with this thruster and this may change the data of 
Figure 52. Curves similar to those of Figure 52 have been obtained at 
higher discharge loss levels and this suggests that, while the entire 
discharge loss curve will be shifted downward and the entire beam current 
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curve will be shifted upward as a result of discharge chamber improvements, 
the curves of Figure 52 do indicate the extent of the performance degradtion 
associated with biasing the screen grid. Additional tests which will focus 
on improving this thruster and verifying this statement are planned during 
the next report period. 

Additional Radial Field Thruster Observations 

During the course of the tests mentioned in the preceding section 
several features of the radial field thruster were observed which warrent 
discussion. First is the tendency for electrons ejected into the dis- 
charge plasma along particular field lines to tend to remain attached to 
those lines. This was evidenced for example in the configuration of 
Figure 50 as a series of luminous spokes of nearly constant diameter ex- 
tending from each hole to the outer wall of the thruster. These spokes, 
which were most visible at high magnetic field intensities, suggest high 
energy electrons interact extensively along the field lines on which they 
are injected into the chamber before they migrate toward the anode. The 
fact that the spokes are of essentially constant diameter suggests that the 
cyclotron radius is small compared to the diameter of the holes, because 
the 1/r variation of the magnetic field would otherwise cause the diameter 
of luminous spoke to increase noticeably near the outer radius of the 
thruster. Numerical calculations of electron trajectories based on actual 
magnetic field measurements in the thruster showed that a spoke of nearly 
constant diameter was predicted as suggested by the filings map of Figure 50. 

A discharge chamber having a configuration like the one shown in 
Figure 53a was made and fitted with a tungsten filament cathode which was 
a circular loop of constant radius located at the region of electron emission 
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identified on the figure. The loop cathode was supported in four places 
by screws affixed to the upstream end of the chamber. A screen at the 
downstream end served as the anode. The device was operated in a bell 
jar and a luminous sheet appeared which was azimuthal ly quite uniform and 
had a cross section similar to that shown by broken lines in Figure 53a. 

The thickness varied with magnetic field strength (as controlled by magnet 
current) but for 10 gauss radial at the outer pole piece it appeared to 
be about 1 cm thick. When two segments of the cathode were removed and 
the thruster operated under the same conditions, the luminous sheet was 
interrupted as shown in Figure 53b. 

A second observation made on the radial field thruster was that it 
was very difficult to establish the- desired impedance between the main and 
cathode discharges with the center pole piece configuration of Figure 50. 

Holes ranging in size from 0.3 cm dia to 0.6 cm dia were used in the pole 
piece and it was found that with hole diameters less than about 0.48 cm 
the main discharge could not be started. With this hole diameter and 
larger, however, the voltage drop necessary for proper thruster operation 
could not be sustained between the main and cathode discharge regions. 

Various numbers of the holes were also closed off in an attempt to increase 
the impedance between the main and cathode discharge regions. Making these 
changes in addition to adjusting the main/cathode flow rate split and up- 
stream anode position (Fig. 50} made very little difference in the performance 
curve knee, which tended to lie around 500 eV/ion at 50% utilization. 

It was determined that the establishment of the desired impedance be- 
tween the main and cathode discharge regions necessitated the incorporation 
of a magnetic field across the apertures. To accomplish this the center 
pole piece configuration of Figure 54 was employed. With this arrangement 
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a magetic field was produced across the baffle aperture in such a way that 
it impeded electron flow from the cathode .discharge region and thereby 
facilitated the maintenance of the desired voltage- drop. Operation was 
attempted over a range of aperture widths and magnet current conditions 
and the thruster performance improved substantially over that observed 
with the configuration of Figure 50. 

In the first design of the type suggested in Figure 54, the pole piece 
tip was supported by three 0.13 cm dia stainless steel (non-magnetic) wires. 
It was observed however that when the aperture width was sufficient to 
achieve the desired aperture impedance the magnetic flux passing from the 
center pole piece tip to the outer pole piece was reduced significantly 
because of the reluctance of the gap. This appeared to allow primary elec- 
tron access to the screen grid when it was biased toward anode potential. 

To correct this, the three stainless steel wires were replaced by three 
ferromagnetic steel wires which saturated during operation so the magnetic 
field required to assure the proper aperture impedance still existed. Of 
several configurations of the type shown in Figure 54 which were tested, the 
one having an aperture width of 0.86 cm and utilizing 0.16 cm dia mild steel 
wire to support the tip gave the best thruster performance (that shown in 
Figure 51). 

Further improvements in performance are considered likely in this 
thruster, through such changes as moving the cathode downstream into the 
vicinity of the aperture and possibly introducing multiple slot apertures 
over the length of the center pole piece. These should reduce losses 
associated with the cathode discharge and improve the uniformity of the 
primary electron distribution over the main discharge region respectively. 
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Thruster Stability 

When this thruster was tested with the .holes serving as the aperture 
(Figure 50) it was observed to operate fairly stably at either low magnetic 
fields where the discharge voltage was too low for good performance, or at 
high magnetic fields when the discharge voltage was excessive. During the 
first tests conducted with the slot aperture (Figure 54) and non-magnetic 
support wires it behaved similarly. . In these instances the discharge tended 
to become more stable as the screen grid was biased closer to the anode 
potential. The stability did not seem to be affected by which magnet 
(radial or axial) was being used to produce the field. In fact both sta- 
bility and performance appeared to be the same when the magnetic field level 
in the discharge chamber reached given values, independent of whether the 
field was produced by the axial magnets, the radial magnets or some combination 
of the two. 

When the slot aperture with ferromagnetic support wires was used, the 
thruster operated quite stably over the full range of magnet currents. Sta- 
•bility was observed however to be enhanced through the use of the axial 
magnets in preference to the radial ones. Biasing the. screen grid toward 
anode potential did not appear to enhance stability in this case. 



THE DIFFUSION AND COLLECTION OF ELECTRONS IN ION THRUSTERS 

Glen Longhurst 

The stream of high energy electrons, which conveys ionization enery 

into the main discharge of electron bombardment thrusters, is accelerated 

into this discharge by the potential difference which exists between the 

cathode region and the main plasma. These primary electrons are generally 

assumed to be monoenergetic and to be confined to a region bounded by the 
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surface of revolution of the critical field line. This line is defined 
as the one passing from the outer edge of the cathode pole piece to the 
anode pole piece. Generally one would expect primary electrons to be 
bound to lines which pass through or very near the aperture where these 
electrons enter the main discharge. For hollow cathode thrusters this 
aperture is between the baffle and the cathode pole piece wall, and it has 
an annular shape like the dotted surface shown in Figure 55a. From this 
figure it can be seen that true .primary electrons (those which have had no 
energy altering collisions) should be confined to the region between the 
surfaces of revolution of the magnetic field lines intercepting the baffle 
aperture at its inner and outer edge. In the model being suggested here 
the primary electrons move back and forth along field lines and are re- 
flected at the baffle aperture and downstream end of the discharge chamber 
by the electric fields that exist at these locations. 

Eventually true primary electrons have collisions, their energies are 
altered, and they can move off of the field lines on which they were located 
initially. Although the electrons may radially move in or out, the net 
drift is radially outward into the "collisional migration region" of 
Figure 55a, which again has its shape determined by the shape of magnetic 
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field lines. This is the region through which the electrons must pass to 
reach the "electron collection layer" from which' they can be taken to the 
anode. The outer surface of this layer is the surface of revolution of the 
magnetic field line that intercepts bhe anode. The electron collection at 
the anode is of course necessary to allow a continuous stream of primary 
electrons to feed energy into the plasma while preserving plasma neutrality. 
The electrons, which pass through the "collisional migration region, con- 
tinue to move readily along field lines while being reflected at end sur- 
faces and they are assumed able to move across field lines only , as a result 
of collisions. With each collision the electrons lose energy, and so one 
would expect to see the primary electron energy distribution function change 
from a monoenergetic one in the "true primary electron region" to a lower 
average energy and move widely spread distribution in the "electron collec- 
tion layer." In this region the electrons probably have a distribution 
function that is essentially Maxwellian. It is important to note that this 
proposed model, the restriction to electron migration is across field lines. 
The processes by which these lower energy electrons migrate to the anode 
and are subsequently removed from the plasma are of interest because they 
govern the discharge current which can be supplied and the rate, at which 
energy supplied by the discharge current is carried to the anode by the 
electrons. 

Plasma Properties along Magnetic Field Lines 

If one can assume that electrons undergo essentially unrestricted 
motion along magnetic field lines compared to that across them, then 
analysis of the electron migration problem is greatly simplified. The 
assumption appears valid because the collision frequency is generally 
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small compared with the cyclotron frequency of electron motion around the 
magnetic field lines. As a result electrons traverse the length of the 
magnetic field line to which they are bound several times before having 
a collision. Further it is observed that where magnetic field lines through 
the plasma intersect only surfaces which are at cathode potential, electrons 
bound to those field lines and arriving at the plasma boundary, will be re- 
flected back into the plasma. This is the case for all of the regions iden- 
tified in Figure 55a except the electron collection layer. There, electrons 
arriving at the anode potential plasma boundary with momentum normal to the 
boundary surface great enough to overcome the sheath potential will be 
collected by the anode and removed from the electron population. 

The relative freedom of electron motion along magnetic field lines 
implies that intensive electron properties should be uniform along field 
lines. To check the validity of the assumption, lines of constant Max- 
wellian electron temperature were plotted and compared with the associated 
magnetic field configuration for three different thrusters. These were 
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the axial and divergent (SERT II) field configurations tested by Knauer 
and the cusped magnetic field (CMF) thruster tested by Beattie. These 
comparisons are shown in Figure 56 where constant Maxwellian electron 
temperature lines are solid while magnetic field lines of force are broken. 
Although the configurations are not congruent, qualitative agreement is 
quite good considering errors inherent to Langmuir probe positioning and 
trace analysis and the linear interpolation of the resulting plasma 
property data. 

Figure 57 shows a similar comparison for the Maxwellian electron 
density. Here constant density lines (solid) are compared with the same 
magnetic field lines (broken). Again the broken and solid lines show 
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reasonable congruency when one considers the fact that a decrease of 
density in regions of diverging magnetic field is expected. This is ex- 
pected because the electrons tend to be bound -to field lines, -and as the 
field lines separate, the density of electrons should diminish. 

Because the intensive electron properties tend to be constant along 
magnetic field lines, it is possible to simplify two-diminsional geometries 
such as the one of Figure 55a to the one-dimensional configuration shown 
in Figure 55b. In the simplified model the magnetic field is assumed to 
be uniform and normal to the plane of the figure. The virtual anode sur- 
face shown in Figure 55a is the surface of revolution of the innermost 
magnetic field line which intercepts the anode. It becomes a circle in 
the one-dimensional model (Figure 55b). The greater mobility of electrons 
along field lines than across them implies that only a few electrons will 
be diverted by collisions to orbit field lines outside the virtual anode 
surface before they are collected by the anode. That this is so is 
supported by the low plasma densities measured outside the virtual anode 
surface and by the observation that nearly all of the anode current is 
drawn to that region of the actual anode where it is closest to the virtual 
anode surface (downstream end in the case of Figure 55a). It is noted in 
passing that thrusters in which the primary electron region overlaps the 
electron collection layer exhibit poor performance due to the loss of 
primary electrons before they can ionize neutral propellant atoms. 

Electron Diffusion through the Collisional Migration Region 

The electron current density J to the electron collection layer 

0 
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should be describable by classical diffusion theory. 
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where e is the electron charge, the electron density, D the diffusivity, 
y the electron mobility, and E the local electric field. Assuming the 
plasma potential gradients are small inside the electron collection layer, 
the steady-state electric field E is negligible. Neglecting the second 
term on the right. Equation 12 becomes for the simplified model of Figure 54b 


e = -D 
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(13) 


Here I„ is the net electron current drawn by the anode. A' is the area 
" A 

of the virtual anode surface, is the ion current density to the actual 
anode as determined from the Bohm criterion for a stable sheath, and e is 
again the electronic charge. Using measured density gradients (dn^/dr) 
across the "collisional migration regions" of the three thrusters considered 
in Figures 56 and 57, values of D required to satisfy Equation 13 were com- 
puted, These are listed as measured diffusivities in Table III. 

Classical theorgy gives the electron diffusivity (D|^) across the weak 
magnetic field of a laminar plasma as 


T m 
_ mx e 

e~W 


(14) 


where T^^^ is the electron temperature (eV), m^ is the electron mass, e is 
the electronic charge, and B is the magnetic field intensity. The electron 
collision frequency v is approximated by 


v=(na +n a +n, 
' e e 0 0 H 


a+) U 


( 15 ) 


where n^, n , n. are the densities of electrons, neutrals and ions, 
6 0 ’^ 
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respectively and where the corresponding a's are total collision cross 
sections with those species for electrons with an average speed U . It 


is assumed that- 


i 



( 16 ) 


Table III. Ion Thruster Diffusivities 



Using measured values of plasma properties near the electron collection 
layer in the above equations, for each of the thrusters under con- 
sideration was calculated. These values are listed as laminar diffusivi- 
ties in Table III. They are observed to be lower than the "measured" 

values by about two orders of magnitude. 

> 

It has been observed that the laminar diffusion model does not 
accurately account for electron current densities observed in plasmas 
similar to those under consideration here.^^’ Bohm^^ has 

suggested that the equation 


( 17 ) 
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which was empirically obtained from a variety of data be used. The Bohm 
diffusivities were calculated for the three thrusters considered, and 
values appear in Table III. The agreement with "measured" values is ob- 
served to be considerably better than for laminar diffusivity calculations. 
According to Bohm, Equation 17 is only accurate to within a factor of 
2 or 3 and this is indeed the sort of agreement between "measured" and 
"Bohm" values in Table III. 

Departure from "laminar" diffusion to what has been called "turbulent" 
45 40 42 43 

diffusion by Bohm and others ’ ’ is generally attributed to plasma 

instabilities. Analyses of such instabilities lead to several different 
modes of plasma oscillations depending on the nature of the plasma and 
associated magnetic fields. One such mode predicts an azimuthal component 
of an oscillating electric field. This component of the electric field 
reacts with the generally axial magnetic field to produce radial drift. 

The modes actually present in thruster plasmas remain to be verified before 
a conclusion can be drawn regarding the mechanism of electron diffusion 
with the plasma. 

Electron Collection at the Anode 

When electrons reach the electron collection layer of Figure 55a they 
become subject to the collection process. The one-dimensional model of 
Figure 55b is again assumed to apply. For mathematical simplicity the 
distance variable r is measured from the virtual anode surface as shown. 

In this single dimension the helical motion of electrons around 
magnetic field lines is observable only as a sinusoidal oscillation in 
time, as suggested by Figure 58. The period of this oscillation would be 

2iTm 

e 

eB 


X = 


(18) 
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_J L T, = ^Trmp 

nr eB 



ONE-DIMENSIONAL ELECTRON MIGRATION TO ANODE 


FIGURE 58 
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Again is the electron mass, e is its charge, and B is the magnetic 
field intensity. The amplitude of the oscillation, which represents its 
maximum excursion from its average location r, is the Larmor radius 


m U Sin 0 
_ e 

“c ea 


(19) 


where U is the electron's speed and e is the angle between the velocity 
vector and the local magnetic field vector. This angle is constant over 
the interval of time between collisions and is not altered by specular 
reflections from cathode potential surfaces. The average location r is 
identically the location of the magnetic field line which is the electron's 
guiding center during some time interval At. As an' electron has collisions, 
the location of its guiding center changes randomly and discontinuously as 
suggested by Figure 58, By defining a differential volume element per unit 
length to be 

d (^) = 2ti (R-r) dr (20) 

where R is the radius of the virtual anode surface, electrons may be said 
to become residents of d(-V-/L) when their average locations r fall between 
r and r + dr. 

If as a result of collisions an electron becomes bound to a field 
line which is at most one Larmor radius away from the virtual anode sur- 
face, then under the assumption that its collision frequency is much less 
than its cyclotron frequency, it will travel along that field line until 
it comes near the actual anode. If it has sufficient momentum to overcome 
the sheath potential which exists there, it will be collected by the anode. 
As a result of this collection process one would expect the electron energy 
distribution function measured within the electron collection layer to be 
deficient in electrons with sufficiently high energies to be collected. 
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The collection process can be quantified if one can determine the 
probability of electron collection in terms of the electron's average 
location T and its kinetic energy e. Let the potential difference between 
the plasma near the anode and the anode itself be V^. It will be assumed 
here that the plasma is more positive than the anode making positive 
(this is generally the case in ion thrusters). is usually 1-5 V. Because 
the Debye length for thruster plasmas is so short, usually of order 10"^ m, 
the sheath thickness may be neglected. Anode curvature may also be neg- 
lected so that the critical distance from the anode r^ to a magnetic field 
line from which an electron with Larmor radius can be collected is shown 
in Figure 59. The radial component of electron momentum when it encounters 
the sheath must be sufficient to carry it across the sheath. This criterion 
for collection then becomes 

■ > ■ m 

“ (U Sin 9 Sin a)^ > eV^ , (21) 

where the angle a as shown in Figure-58-is given by 
a = Cos"^ ^ 

Substituting Equation 22 into Equation 21 gives 


m L|2 

-i — Sin^e 



> e . 


(23) 


m 

6 

The greatest distance r^,„ from which an electron having energy e = — 

niaX ^ 

can be collected corresponds to the case when with speed U the electron 
is moving at right angles to the field vector (sin 0=1), and from 
Eqn. 23 it is given by; 


I 






ANODE 

SURFACE 



ELECTRON COLLECTION AT ANODE 
FIGURE 59 
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If one uses the maximum velocity at which 'significant numbers of electrons 
are present, then this equation defines the thickness of the electron 
collection layer. 

m 

If one now considers electrons having an energy e = -| — which are 

at. a location r (less than r ) then Equations 19 and 23 combine to 

max 

reveal that only those electrons having an angle e between their velocity 
vectors and the field vector given by the following will be collected: 




(25) 


Assuming an isotropic distribution of velocity vectors at all values of r, 
the fraction of electrons with energy s which will be collected from 
location r is called f(er) and it is given by 


f(e,r) 



2ir Sin e de 
4tt 


* 

Cos 0 , 


(26) 


where e represents the minimum angle between the velocity vector and the 

magnetic field vector for which collection is possible. From Equation 25 
* 

0 is given by ± 


■k 

e 


Sin-^ 


(erB)‘ 

2em. 


e V, 


+ 



(27) 


Using. Equation 26 and letting 



(erB)^ 
2 m 


+ e V 

s 


e 

be the minimum kinetic energy for collection at average location r 
magnetic field B with sheath potential V^, one obtains 


.( 28 ) 


in 
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f(s,r) 



* 

e > e ■ 


* 

e < e 


(29) 


In the event that is negative, all electrons arriving at the- -anode 

sheath would be attracted to the anode and hence would be .collected. In 

that case Figure 59 would" show a = 0 and r = p . This causes the term 

c c 

e in Equation 28 to vanish. 

The application of Equation 29 to the electron collection process 
can be best illustrated by an example. Suppose that in a particular 
thruster, electrons in the vicinity of the electron collection region were 
described by the Maxwellian distribution function 


g(e) = 


vV eT 


mx 


eT, 


m 


exp 




(30) 


If the temperature of these electrons was 5 eV then the distribution of 
electrons would be as shown by the solid line^of Figure 60. At locations 
close to the virtual anode surface, however, a fraction f(e,r) of the 
electrons would be removed from the distribution function so the distri- 
bution function of the remaining electrons would be given by 


g(ei.r) = g(e) [l-f(e,r)] . (31) 

Assuming now that = 0 and B = 13 gauss, e (r), f(e,r), g(e) and g(e,r) 
can be calculated from Equations 28, 29, 30, and 31 respectively for various 
values of r. The results of doing this are shown in Figure 60. Infinitely 
far in from the anode (r = «>) the distribution function is Maxwellian. At 
a distance 1 cm from the anode, the distribution function is Maxwellian for 
energies less than about 14 eV and beyond this point some of the electrons 
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ELECTRON DISTRIBUTION FUNCTION NEAR ANODE 


FIGURE 60 
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in the high energy tail have been removed as suggested by the dashed line. 

As one gets closer to the anode (r = 0.75 and 0.5 cm) more of the elec- 

* 

trons are removed and the cutoff energy e extends to lower values. 



SIZE/SPECIFIC-IMPULSE SCALING RELATIONSHIPS 
H. R. Kaufman 

In- view of the fact that ion thruster costs tend to scale with the 
number of components rather than their size, it is probably more cost 
effective to use one large thruster rather than several smaller ones. 
Additional thruster cost savings might be achieved by going to very high 
specific impulses to increase thethrust density and therefore the thrust 
■of a unit of specified size. The extent to which thruster sizes and 
specific impulses can be increased under the constraints of acceptably 
long lifetime and good thruster performance is the subject of this analysis. 

Analysis 

Assume Maxwellian electron temperature, primary electron energy, and 
the primary-to-Maxwell ian electron density ratio are all constants. These 
assumptions should.be reasonable for operation near the "knee" of the dis- 
charge-loss/propellant-utilization performance curve for thrusters with 
diameters greater than about 30 cm. Since there is a tendency for electron 
temperature to decrease with thruster diameter this assumption is con- 
servative. 

28 

Using the above assumptions and the doubly charged ion model , one 
obtains 

where n^^ is the density of doubly charged ions, n^ is the total electron 
density, and Vp and Ap are the volume and surface area of the primary 
electron region. Inasmuch as n.^.^ « n^ is required to assure a long 
screen grid lifetime, n^ = n^ and the above becomes 
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n. ,/n. “ n V„/A^ 
++' + e p p 


(33) 


From Equation 33, it can be concluded that the value of n^V /A„ should 

e p p 

not be much greater (if any) for larger thrusters than it is for the 30 cm 
thruster — otherwise discharge-chamber erosion would become more serious. 
Keeping this parameter constant, though will limit the propellant utiliza- 
tion which is given by the equation 


J. /J 
b 0 


\ ■ T ¥ J JT 


(34) 


where is the ion beam current density, and is the loss rate of 
un-ionized neutrals (same units as Jj^). Equation 34 clearly shows that 


db/Jp should be as high as possible to maximize propellant utilization. 
For a given propellant. 


(35) 


46 

where n^ is the neutral density in the discharge chamber. Under the 
constraint of constant electron temperature, energy and primary- to- 
Maxwellian density ratio, the ion-beam current density becomes 


J, a n n V /A 
b e 0 p p 


(36) 


The ratio with substitution of Equations 35 and 36 is, then. 


J, /J 04 n V /A 
b o e p p 


(37) 


But Equation 37 uses the same grouping n^Vp/A^ that determines the 
ratio n^^/n^. Limitation of the doubly-charged ion ratio thus leads to 
a similar limitation in maximum propellant utilization. 
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, The maximum value- of n V /A can also be used to infer a maximum 

e p p 

g 

depth for the- discharge chamber. From the Bohm velocity criterion the 
ion-beam, current density can be related to the discharge-chamber electron 
.density, 

J, a n (38) 

b e 

Combining Equations 36 and 38, and solving for Vp/A^ yields 

Vp/Ap-l/Ho . - (39) 

Since n V /A is limited to a maximum value however Equation 36 becomes 
e p p 

■>5 ‘ "o ’ 

and Equation 39 can therefore be rewritten as 

V /A - i/j . (41) 

p p b 

For -a large, flat discharge chamber, the depth is about twice V /A . 

r r 

This depth,:, then will vary inversely as ion-beam current density. 

The energy cost of making ions in the discharge tends toward a con- 
stant of 50-60 eV/ion.^^ The higher loss per beam ion observed in actual 
thrusters. is explainable in terms of ion-chamber geometry. For large 
flat chambers, between 1/3 to 1/2 of the ions produced should leave the 
chamber via the ion beam. Discharge losses should therefore be in the 
100-180 eV/ion range for large, flat chambers. 

The maxi mum. useful thruster diameter can be determined from breakdown 
voltages and acceptable span-to-gap ratios. For the 1-10 mm range of grid 
gaps with mercury, the maximum allowable electric field is roughly 

c nn 

2 X 10°V/m.^ For a reasonable net-to-total voltage ratio R of 0.7, the 
net accelerating voltage Vpet would be about 1400 V per mm of accelerator 
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system, gap. Assuming dished molybdenum grids the beam diameter should not 

47 

exceed 300-600 times the gap. For singly ionized mercury, the specific 
impulse is approximately 


sp net 


(42) 


Substitution of 1.4 x 10 v/m times the gap for gives 


= 100(1.4 X lo'^a . 
sp g' 


(43) 


Solving for i yields 
y 


£ = 7 X 10“^'' I 

g sp 


(44) 


Using 300-600 times this gap as the beam diameter one obtains^ 


D. = 2 - 4 X lO'^I ^ 
b sp 


(45) 


The permissible beam diameter, from structural and breakdown considera- 
tions, thus increases as the square of the specific impulse. If the 
span-to-gap ratio of the accelerator system is held constant, and R is 
not decreased, an increase in diameter over that given by Equation 45 
will result in no increase in thrust. 

The charge-exchange erosion will also be reasonable for higher 
specific impulses. This can be shown .from Child's law. 




(46) 


where V is the total accelerating voltage. For a constant electric field 
strength E, between- the grids this becomes 

cc E^/V^ « 1/v^ (47) 
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The sputtering yield Y varies as where n is a number between 0 and 
1 over the energy range of interest. The rate of accelerator material 
removed per unit area m,^ is given by 

aC 


m J. n Y £ . 

ac bog 


(48) 


With the substitutions of for Y and for n^ this becomes 


"ac “ ’‘g • 


(49) 


Noting that will tend towards a constant for an optimum design, 

using Equation 47 to substitute for and recognizing a constant 

electric field implies I V, Equation 44 becomes 

y 

lilac " 


Charge exchange erosion per unit beam area will therefore increase as 
specific impulse is increased -but compensation for the increased erosion 
can be achieved by making the accel gHd thickness proportional to 2,^. 

Designing for a constant electric field and a constant net-to-total 
voltage ratio also results in' constant thrust per unit beam area. This 
can be shown from 

T/''“Jb'sp ■ ’ <51) 

Substitution of Equations 42 and 47 and holding R constant gives 


T/A = constant . 


Conclusions 

The conclusions which can be drawn from a study such as this depend 
on the level of technology used in the analysis. The level assumed here 
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has been in general, that of the current ’30-cm mercury thruster. Future 

technology developments could obviously affect the validity of the con- 
clusions given. 

1. To keep the doubly ionized fraction within reasonable limits for dis- 
charge chamber erosion, the maximum propellant utilization will be , 
limited. Based on present 30 cm data, this limit will be 0.90-0.95. 

2. The limits on doubly ionized fraction and propellant utilization 
also imply a definite depth for the discharge chamber. This depth is 
small enough that the discharge chamber for a large, high specific 
impulse thruster will have L/D « 1. To maintain a uniform beam with 
such a chamber (without undue complexity) would require a low magnetic 
field strength over most of the volume. A multipole magnetic field 

is one configuration that would meet this low field strength requirement. 
From 30 cm data, the chamber depth should be 5 to 10 cm at 3000 seconds 
specific impulse and from Equations 41, 42, and 47 it should increase 
in proportion to the I^p at higher values. 

3. The fraction of ions produced that escapes- into the ion beam should 
increase with thruster size and this should result in a lowering of 
the discharge losses (100-180 eV/ion). 

i 

4. The maximum useful thruster diameter increases rapidly with specific 
impulse. From Equation 45 one obtains the following: 

L^p(sec. ) 3000 5000 7000 10000 

Dj^(cm) 18-36 50-100 100-200 200-400 

5. The thrust per unit beam area should be approximately constant over 
a range of specific impulse. A factor of perhaps two improvement 
might be expected over 30 cm data— from a more uniform beam profile. 

6. Under the imposed constraints of a -fixed doubly- to-singly charged 
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ion fraction in the discharge chamber and a constant electric field 
between the grids, the- discharge-chamber erosion rate will decrease 
with increasing specific impulse. The accelerator erosion rate will 
increase with specific impulse, but can' be offset by making grid 
thickness proportional to accelerator-system gap. Because of these 
effects, there may be a net advantage to increasing the doubly ionized 
fraction (and hence propellant util izati on) slightly at high specific 
impulses. 

7. Multiple discharge-chamber cathodes will probably be required for 
thrusters much larger than 30 cm in diameter. The problems of 
parallel hollow cathodes indicate that it would be simpler to use a 
separate discharge supply and control loop for each cathode. 



MOLECULAR NITROGEN ION DENSITIES IN 


MERCURY DISCHARGE CHAMBERS 

During the period covered by this report anomalies in the erosion 
rates of cathode potential surfaces in mercury ion thrusters were dis- 
covered in tests conducted at two different laboratories. An experiment 
was performed that suggested these anomalies might be caused by variations 
in the background partial pressure of nitrogen in the test facilities. It 
was postulated that this occurs because some of the nitrogen molecules 
entering the discharge chamber are ionized. Subsequently they strike the 
truster walls and because they are in an ionized, reactive state they form 
a sputter-resistant material with the wall. If this process does occur it 
might provide an explanation for the experimentally observed increase in 
putter resistance of thruster materials with nitrogen partial pressure. 

One step associated with the study of this phenomenon was the de- 
termination of the density of molecular nitrogen ions in the discharge 

chamber of a thruster. To calculate this quantity the model deve-loped 

28 

for the calculation of doubly charged mercury ions was adapted to the 
problem as described below. Nitrogen ion densities were then calculated 
for the operating conditions of the 30 cm diameter thruster in which the 
anomolous erosion rates had been observed. 



Figure 61 is a schematic of a vacuum tank containing nitrogen at a 


partial pressure and a density n^ together with a discharge chamber. 

If the vacuum tank surfaces are at an average temperature T^ (°K), then the 
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NITROGEN FLOW SCHEMATIC 


FIGURE 61 
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nitrogen density in the tank will be given by 


"o ' kf: 


(52) 


where k is Boltzmann's constant. Some of these neutral nitrogen molecules 
flow into the discharge chamber at a rate n^ determined by the equation 
for free molecular flow through a sharp-edged orifice. 


V A. 
0 0 A 


(53) 


Where A^ is the open area of the grid system to neutral flow (generally the 
open area of the accelerator grid) and v^ is the thermal velocity of the 
nitrogen. 


V 

0 



(54) 


The term m^ in this equation is the mass of a nitrogen molecule. 

The nitrogen flow into the thruster as given by Equation 53 must 
be equal to the nitrogen outflow in the steady state. This outflow is 
the sum of the neutral outflow rate and the ionic outflow rate and is given 
by 

n V A„ 

n = v^ A 3 (55) 


where n is the neutral nitrogen molecular density inside the discharge 
chamber, n^ is the corresponding molecular nitrogen ion density and A 3 is 
the open area of the grids to ions (generally close to the screen grid 
open area). The term v in this equation is the velocity of the neutral 
nitrogen molecules' within the discharge chamber and is given by 


/8kT 


w 


irm. 


V = 


(56) 
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where T, , is the discharge chamber wall temperature (°K). The term v, in 
Equation 56 is the nitrogen ion velocity as determined by the Bohm 
equation 48 

V 

As sugge'sted by this equation the electrons effecting the ionizations 

in the discharge chamber are assumed to consist of a Maxwellian group 

(temperature T (eV) and density n ) and a monoenergetic or primary group 
mx mx 

(energy ?p^(eV) and density Hp^). The term "e" in Equation 57 is the 
electron charge. 

In order to calculate the ion density for use in Equation 55 it was 
necessary to model the ionization process occurring in the discharge chamber. 
It was therefore assumed that only one reaction was important in the pro- 
duction of molecular nitrogen ions, namely 

N 2 + e~ ^ + 2e” . 

49 

Cross sections for this reaction were obtained from the literature and 
other possible reactions involving intermediate dissociation and/or elec- 
tronic exitation reactions were neglected. These other reactions could 
be significant but their incorporation into the analysis would require a 
thorough literature survey which was not considered to be warrented here. 

The process defined by the above equation should produce results accurate 
to within t 100%. 

If one equates the production and loss rates of nitrogen molecules in 

28 

the manner described by Peters the following equation for the ratio of 
neutral-to-sing,ly ionized nntrogen molecules in the discharge chamber is 
obtained. 
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The previously undefined terms appearing in this equation are: 

V /A - vol ume-to-surface area ratio for the primary electron region 

r r 

of the discharge chamber. 

- the uniformity factor describing the ion flux averaged at 
the surface of the primary electron region to the flux averaged 
over the volume of this region. 

4 - 

Po (6py.) " The cross section - electron velocity product pertaining to 

primary electrons at the energy This parameter is plotted 

in Figure 62 for the reaction being considered. 

Q (T ) - The cross section - electron velocity product integrated 

U lilX 

over energy for Maxwellian electrons at the temperature T^^. 
This parameter is also plotted in Figure 62 for the reaction 
being considered. 

Combining Equations 52 through 57 one obtains 



or rearranging to solve for th.e nitrogen ion density 



Using the geometrical and plasma conditions pertaining to the particu- 
lar thruster being considered Equation 58 can first be solved and then 
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FIGURE 62 
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Equation 60 can be solved for the nitrogen molecular ion density. 

Application to the 30 cm Thruster 

For the 30 cm thruster with EM optics operating at a 2 amp beam current 

' oo 

the following properties apply: 

V-/A = 0.025 m 
F+ =1.5 

n^^ = 3.28 X 10^® m“^ 
pr 

n_ = 13.1 X 10^^ m"^ 
mx 

V = ^'-2 

n /n =0.25 
pr mx 

^ ^ 0.676 ^ -I .f. 

A^ 0.433 ^ ^ 

+ -13 3 

For these conditions Figure 62 gives = 0.25 x 10 m /sec and 


+ - ! b ^ 

= 0.4 X 10 m /sec. If the temperature of the tank liner is assumed 


^-15 3 


to be 


Tq = 100°K , 


and the temperature of the thruster wall is assumed to be 


then Equation 58 gives 


T = 500°K , 
w ’ 


ji 


123. 


n, = 1.98 X 10^® P 
+ 0 


Equation 60 then gives 
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2 + -3 

where in nt/m yields n in m . Substituting values of P^ in the range 
of interest the results shown in Table IV are obtained. 

Table IV 

Molecular Nitrogen Ion Density in the 30 cm Diameter EM Thruster 


P (torr) 

n^ (cm“^ 


2.6 X 

10 

10"^ 

2.6 X 

10' 

10-s 

2.6 X 

10 

10"^ 

2.6 X 

10 
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